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2.  EXECUTIVE  SUMMARY 

Toward  replacing  current  neutron  detection  technologies,  solid-state  designs  offer  the  promise 
of  smaller/lighter  detectors  with  lower  power  needs  and  monetary  costs.  Direct-conversion 
detectors,  in  contrast  to  conversion-layer  detectors,  are  capable  of  achieving  100%  thermal  neutron 
detection  efficiencies  with  thermal  insensitivity  and  lower-complexity  heterostructure  designs. 
However,  as  these  are  based  on  immature  semiconductor  technologies,  they  remain  in  the 
fundamental  research  stage.  A  number  of  authors  have  reported  the  successful  fabrication  of  direct- 
conversion  neutron  detector  heterostructures  from  amorphous  hydrogenated  boron  carbide  (a- 
ByC:Hv)  grown  on  n-Si  by  plasma-enhanced  chemical  vapor  deposition  (PECVD)  from  carborane 
precursors.  These  have  been  limited  to  devices  ~1  pm  thick  with  detection  efficiencies  of  <1%.  It 
has  been  found  that  the  experimental-to-theoretical  efficiency  drops  as  the  thickness  of  the  a- 
B.yC:Hv  layer  increases,  which  highlights  the  tradeoff  between  higher  neutron  capture  efficiency 
from  thicker  films  and  decreased  charge  collection  efficiency  owing  to  the  poor/unoptimized 
charge  transport  properties  of  this  material.  The  objective  of  this  effort  has  been  to  rigorously 
investigate  the  charge  transport  properties  of  a-BvC:Hv  films  to  understand  how  both  interfacial 
and  bulk  transport  contribute  to  detector  he tero structure  performance,  while  seeking  to  optimize 
these  properties. 

Over  a  hundred  a-BxC:Hy  films  were  grown  with  varying  PECVD  conditions.  Films  were 
deposited  with  an  extremely  wide  range  of  properties,  including  growth  rate  from  1-1 12  nm/min, 
B/C  ratio  from  2. 8-5.0,  atomic  concentration  hydrogen  (at.%  H)  from  10-45%,  density  from  0.9- 
2.3  g/cm3,  hardness  from  1-25  GPa,  Young’s  modulus  from  10-340  GPa,  Tauc  optical  band  gap 
from  E7-4.0  eV,  Tauc  slope  parameter  {BU2)  from  -200-2200  cm_1/2,eV_1/2,  Urbach  energy  from 
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0. 1-0.8  eV,  dielectric  constant  from  2. 4-7. 6,  and  electrical  resistivity  from  108 — 1015  Q-cm.  Strong 
correlations  were  found  between  chemical  composition  (mainly  atomic  concentration  hydrogen, 
and  to  some  extent  atomic  concentration  carbon)  with  mechanical,  optical,  and  electronic 
properties. 

The  chemical  and  electronic  structure  of  the  a-B,vC:Hv-to-metal  interface  was  studied  using 
photoemission  techniques  for  six  different  metals:  Cr,  Ti,  Al,  Cu,  Ni,  and  Pd.  It  was  found  that  Cu 
and  Ni  form  clean  interfaces,  whereas  the  remaining  metals  form  interfacial  oxides.  Based  on  band 
lineup  arguments,  it  was  concluded  that  Cu  should  form  a  reasonable  ohmic  contact,  but  that  Ni 
or  Pd  may  be  better  suited  as  a  hole  injection/electron  blocking  electrode.  The  current-voltage 
study  of  metal-insulator-metal  heterostructures  with  various  metals/semiconductors  revealed  that 
all  of  the  materials  tested  (including  Cu,  Al,  Pd,  Ti,  and  Si),  regardless  of  work  function  or 
interfacial  oxide  formation,  formed  ohmic  contacts.  This  was  rationalized  in  terms  of  the  very  high 
resistivity  of  a-BrC:Hv  (and  thus  ease  of  attaining  ohmic  injection)  as  well  as  its  amorphous 
structure,  with  localized  states  within  the  gap  capable  of  facilitating  charge  injection. 

To  study  bulk  transport  properties,  both  steady-state  and  transient  (dark-injection)  space- 
charge-limited-current  (SCLC)  techniques  were  applied  to  measure  mobility,  and  steady-state 
photoconductivity  techniques  were  applied  to  measure  mobility-lifetime  product.  Based  on  the 
steady-state  SCLC  measurements  of  a  large  number  of  samples,  it  was  found  that  mobility  in  a- 
B.yQHv  ranged  across  ten  orders  of  magnitude,  from  10-5— 10-14  cnr/V*s,  resistivity  from  108-1015 
Q-cm,  and  charge  carrier  concentration  from  1014-1018  cm-3.  Transient  SCLC  measurements  were 
completed  on  a  subset  of  samples,  and  values  determined  using  the  technique  were  consistently 
higher  than  those  determined  using  the  steady-state  analog.  This  was  rationalized  as  a  consequence 
of  slow  trapping,  which  is  reflected  in  the  longer  time  scale  of  the  steady-state  experiment.  Overall, 
the  transient  experiment  may  be  more  relevant  to  neutron  detection  as  it  is  probing  a  similar 
timescale.  Mobility-lifetime  product  (jut )  measurements  performed  on  low-mobility  samples 
yielded  jit  values  on  the  order  of  10'13-10-14  cm2/V,  in  the  same  range  as  the  mobility,  which 
suggested  a  relatively  long  lifetime  on  the  order  of  1  s.  Mobility-lifetime  product  determinations 
for  higher  mobility  samples  and/or  direct  measurements  of  carrier  lifetime  will  be  needed  to  verify 
the  generality  of  this  observation.  Charge  transport  properties  did  not  show  a  simple  correlation 
with  either  hydrogen  concentration  or  electronic  structure  parameters,  but  rather  showed  two 
distinct  regimes.  This  was  rationalized  as  a  transition  from  a  variable  range  hopping  mechanism 
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for  films  with  low  band  gap  and  wide  band  tails  to  a  multiple  trapping  mechanism  for  films  with 
high  band  gap  and  narrow  band  tails. 

Modeling  of  charge  collection  efficiency  and  neutron  detection  efficiency  for  the  specific  case 
of  a  planar  thin-film  detector  with  low  mobility  was  completed  to  understand  what  charge  transport 
figures  of  merit  were  required  in  a-BAC:Hv  to  achieve  a  target  detection  efficiency.  A  numerical 
method  was  implemented  to  model  the  effects  of  transport  parameters,  mobility,  lifetime,  and 
mobility-lifetime  product,  as  well  as  experimental  parameters,  film  thickness,  electric  field,  and 
integration  time.  It  was  shown  that,  because  in  a  low  mobility  material,  the  transit  time  of  charge 
carriers  may  be  long  relative  to  the  integration  time,  it  is  not  only  / n  that  determined  charge 
collection  efficiency,  but  also  the  individual  n  and  r  parameters.  Further,  because  boron-rich  BAC 
has  such  strong  neutron  absorption  and  very  high  detection  efficiencies  can  be  obtained  with  thin 
films  (10-100  pm),  the  charge  transport  requirements  are  forgiving.  For  a  model  detector  of  10B- 
enriched  B4C  with  a  thickness  of  15  pm  and  excited  carrier  lifetime  of  1  s,  under  an  applied  electric 
field  of  1  MV/cm  and  with  an  integration  time  of  1  ms,  the  mobility  required  to  achieve  the 
maximum  theoretical  detection  efficiency  (-35-40%)  is  5  x  10“ 7  cnr/V-s.  Based  on  the  charge 
transport  measurements  of  a-BAC:Hv  films,  which  demonstrated  that  mobilities  within  this  range 
are  achievable,  it  was  concluded  that — despite  the  poor  transport  properties  of  a-BAC:Hv  in 
general — it  should  be  feasible  to  fabricate  thick  (10-50  pm)  detectors  with  high  efficiencies. 

Overall,  this  effort  has  produced  a  large  body  of  work  investigating  the  charge  transport 
properties  of  a-BAC:Hv  films,  including  the  experimental  techniques  and  analysis  applied  to  extract 
these  results.  This  work  sets  the  stage  for  future  work  to  not  only  demonstrate  that  detection  is 
possible  from  a  boron-rich  direct-conversion  detector,  but  that  neutron  detectors  can  be  made  with 
efficiencies  competitive  with  the  state-of-the-art. 

3.  INTRODUCTION 

3.1.  Amorphous  Hydrogenated  Boron  Carbide  for  Neutron  Detection 

In  the  quest  to  replace  golden  standard  (i.e.,  3He)  and  less  desirable  (e.g.,  10BF3  or  low-flash- 
point  liquid  scintillators)  neutron-sensitive  materials  and  transduction  methods,  1700  US  patents, 
over  3100  active  US  patent  applications,  and  approximately  2700  journal  articles  have  been  issued 
or  published  over  the  last  ten  years.  Of  these  reported  concepts,  three  real-time  detector  class  types 
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have  emerged  in  contention  for  commercial  viability:  (1)  scintillator-based;  (2)  10B-coated  gas 
tube/straw-based;  and,  (3)  semiconductor-based  materials  and  transduction  methods.  Of  the 
scintillator-based  methods,  viable  embodiments  include:  (a)  plastic  scintillators  with  pulse  shape 
discrimination;  (b)  elpasolite  materials;  and  (c)  6LiF:phosphor-based  screen  or  readout  fibers.  In 
all  scintillator  cases,  problems  with  light  readout,  volume,  intrinsic  efficiency,  and/or  production 
cost  have  prevented  these  materials  from  being  accepted  as  industry  standard  replacements. 
Advanced  versions  of  the  10B  gas  tubes/straws  are  serious  contenders  for  large-area 
instrumentation  with  the  benefits  of  position  sensitivity,  but  are  intrinsically  limited  in  their 
efficiency: volume  metric  and  not  of  sufficiently  lower  cost  than  3He.  Of  the  semiconductor-based 
methods  are  indirect-  and  direct-conversion  neutron  detectors.1  Modern  indirect-conversion 
neutron  detectors  with  microstructuring  have  shown  great  promise  with  verified  thermal  neutron 
detection  efficiencies  for  a  single  die  as  high  as  26%2  (or  as  high  as  42%  for  a  double-stacked 
device3,4),  with  controversial  reports  of  efficiencies  up  to  48%. 5,6  However,  processing 
requirements  and  low-band-gap-based  thermal  sensitivity  have  thwarted  wide  acceptance  of  these 
silicon-based  devices.  The  solution  requires  a  reduction  in  the  processing  complexity  of  the 
devices  and/or  reproducing  the  microstructuring  with  higher  band-gap  materials  that  have 
equivalent  electrical  carrier  transport  properties.  Direct-conversion  devices  on  the  other  hand  have 
had  no  commercial  consideration  and  remain  in  the  fundamental  research  category  despite  the 
ability  to  theoretically  produce  such  devices  with  low  complexity  and  with  thermal  insensitivity 
(i.e.,  using  wide-band-gap  materials)  to  achieve  unity  thermal  neutron  detection  efficiency. 

The  challenge  to  realizing  a  direct-conversion  neutron  detection  device  is  that  there  exist  very 
few  semiconductor  materials  that  contain  isotopes  with  adequately  high  neutron-capture  cross 
sections  (3He,  6Li,  10B,  113Cd,  157Gd,  235U/238U),  and  the  few  that  are  known  are  not  well-developed 
from  a  semiconductor  technology  standpoint.  It  is  generally  agreed  that  semiconductors  containing 
6Li  or  10B  are  best-suited  as  solid-state  neutron  detectors  due  to  their  combination  of  high  neutron 
and  low  gamma-ray  cross  sections.2’7  The  class  of  boron-rich  solids  in  particular  boasts  a  number 
of  semiconductor  variants  potentially  viable  for  neutron  detection,  including  BP,8  BN,9  12  B.rC, 13,14 
B12AS2,15  and  Mg2Bi4,16  among  others,  and  various  efforts  have  been  put  forth  to  develop  these 
materials  toward  bulk  crystalline,  composite,  as  well  as  thin-film  solid-state  detectors.  Along  with 
BN,  B.vC  has  been  the  most  highly  studied  class  of  boron-rich  solids  for  neutron  detection.  In 
particular,  BVC  has  shown  promise  as  a  high-resistivity  (>1010  fTcm)  carborane-based  PECVD- 
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grown  a-BvC:Hv  (x  =  3-5)  film  (single  crystalline  BAC  displays  inadequately  low  resistivity  values 
on  the  order  of  10  Q*cm,17  which  implies  problematic  high  leakage  currents  in  a  detection  device, 
as  is  the  case  for  many  other  possible  BVC  variants).  However,  a-BAC:Hv  prototype  detectors  have 
only  been  fabricated  with  very  thin  films  (on  the  order  of  ones  of  microns),  which  implies  that 
charge  pulses  (above  noise)  can  be  obtained  despite  potentially  poor  transport  properties,  but  at 
these  thicknesses  even  ideal  thermal  neutron  detection  efficiencies  are  very  low  (less  than  1  to  a 
few  percent).  Film  thicknesses  in  these  devices  will  need  to  be  on  the  order  of  tens  to  hundreds  of 
microns1218  in  order  to  approach  maximal  efficiencies  for  thermal  neutrons.  A  major  limitation  to 
scaling  up  thickness  are  poor  charge  transport  properties,  which  lead  to  incomplete  electron-hole 
pair  (ehp)  sweepout.  Although  amorphous  BAC:HV  may  be  predicted  to  have  particularly  poor 
transport  properties  due  to  its  disordered  nature,  it  possesses  the  advantage  of  having  a  large 
accessible  configuration  space  of  stable  stoichiometric  and  structural  variants  and  therefore  high 
potential  for  tuning  charge  transport  and  other  material  properties. 

To  ascertain  the  viability  of  a-BAC:Hv  for  thick,  high-intrinsic-efficiency  thermal  neutron 
detectors,  it  must  be  determined  whether  this  material  can  possess  the  right  combination  of  charge 
transport  properties,  including  a  high  resistivity,  to  minimize  leakage  current,  good  electrical 
stability  (i.e.,  a  high  breakdown  voltage  and  therefore  ability  to  withstand  moderate -to-high 
electric  fields),  and  a  sufficiently  high  mean  carrier  drift  length  (proportional  to  mobility  and 
lifetime  via  X  =  /urE,  where  E  is  electric  field),  to  ensure  maximum  charge  collection  efficiency  for 
ehp’s  that  are  generated.1'1  Determining  the  charge  transport  properties  of  a-BAC:Hv  is  complicated 
by  its  high-resistivity,  complex  and  unusual  electronic  structure,  and  disordered  nature.  Many 
traditional  methods  for  measuring  semiconductor  properties  cannot  be  used  in  the  case  of  high- 
resistivity  or  low-mobility  solids.  Analyzing  a  high  resistance  thin  film  pushes  the  limits  of 
experimental  capabilities,  often  requiring  extremely  high  input  impedance  electronics,  additional 
shielding,  and  sophisticated  circuitry  to  surmount  issues  including  very  high  input  resistances, 
extremely  low  current  outputs,  and  complicating  RC  delays.  Further,  even  when  experimental 
measurements  are  feasible,  the  interpretation  of  results  is  complicated  by  the  fact  that  nearly  all 
theory  developed  for  this  purpose  is  based  on  assumptions  (e.g.,  number  and  type  of  defect  levels) 
that  only  hold  true  for  less  complex  materials.  Many  of  the  methods  applied  here  to  study  the 
charge  transport  of  a-BAC:Hv  have  been  developed  and/or  tailored  over  recent  decades  to  study 
first  a-Si:H,  and  later  more  complex  materials  such  as  organic  semiconductors,  and  include  those 
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where  signal  is  actually  enhanced  for  low-mobility  materials  such  as  photoconductivity-based 
methods  (where  low  dark  currents  imply  low  background  noise)  and  time-resolved  drift-mobility- 
based  methods  (where  long  transit  times  lead  to  more  accessible  experimental  timescales).  Further, 
to  accurately  measure  charge  transport  properties  as  well  as  to  fabricate  neutron  detectors  with 
maximal  efficiencies,  establishing  appropriate  electrical  contacts  is  critical.  Depending  on  the  type 
of  device  (e.g.,  p-n  or  schottky  diode)  or  the  type  of  electrical  measurement  (e.g.,  time-of-flight 
drift  mobility  or  Hall  mobility),  either  an  injecting  contact  or  blocking  contact  may  be  required. 
Other  requirements  may  also  include  transparency,  resistance  against  diffusion,  and  environmental 
stability  (resistance  to  thermal  treatment,  oxidation).  Without  having  a  good  understanding  of  the 
electronic  structure  and  charge  transport  at  the  metal -to-a-B.*C:Hv  interface,  one  cannot  be 
confident  in  any  measurements  of  the  bulk. 

Prior  to  this  work,  only  a  few  reports  on  the  charge  transport  properties  of  a-B.vC:Hy  were 
available.  Lee  et  al.  reported  resistivities  on  the  order  of  ~10l()  Dcm  for  similar  a-BrC:Hy  films 
grown  by  PECVD  from  borane  derivatives  in  the  early  1 990 ’s, 20,21  as  part  of  some  of  the  earliest 
work  prompting  interest  in  this  material  for  neutron  detection  applications.  Several  publications 
have  also  reported  low  leakage  currents  (recently  10”5  A/cm2)14  in  reverse-biased  p-n  diode 
geometries.  Although  it  has  been  evident  for  some  time  that  carborane -based  a-BvC:Hv  has  a  very 
high  resistivity,  prior  to  this  work,  no  rigorous  studies  had  been  reported  that  investigated  the  range 
of  charge  transport  properties  possible  in  this  material.  Little  else  is  known  about  the  charge 
transport  properties  of  this  solid,  although  based  on  its  extremely  high  resistivity,  it  can  be  assumed 
that  it  must  exhibit  both  a  low  mobility  and  low  free  carrier  concentration.  a-BAC:Hv  appears  to 
generally  behave  like  a  p-type  semiconductor  based  on  I-V  properties  of  diodes  constructed  from 
this  material,  as  most  evident  in  a  publication  by  Hong  et  al.22  In  a  recent  publication,14  Hong  et 
al.  proposed  a  carrier  density  of  1012  cm”3  and  a  mobility  range  of  10  2-10”4  cm2- V'-s”1  for  a- 
B.YC:Hy  films — values  that  were  not  confirmed  by  measurement,  but  are  not  inconsistent  with 
predictions. 

The  following  report  describes  efforts  to  both  understand  the  charge  transport  properties 
needed  in  an  a-B  vC:Hv  neutron  detector  heterostructure  to  achieve  competitive  efficiencies,  as  well 
as  the  charge  transport  properties  achievable  by  varying  growth  conditions.  Although  previous 
work  on  a-BvC:Hy  has  focused  on  p-n  diode  detectors,  this  work  emphasizes  a  more  simple  planar 
metal-insulator-metal  resistive  detector.  It  is  shown  that  charge  collection  and  detection  efficiency 
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in  a  realistic  experimental  scenario  must  be  modeled  differently  for  thin  film  heterostructures  than 
for  a  typical  detector  based  on  a  thicker  and  higher  mobility  material.  The  techniques  and  analyses 
used  for  investigating  the  electronic  structure  and  charge  transport  of  a-BvC:Hv  are  described,  the 
range  of  properties  measured  are  reported,  and  correlations  between  them  discussed.  This 
information  provides  the  basis  for  tailoring  a-BAC:Hv  films  specifically  for  neutron  detection 
applications  based  on  growth  conditions  and  chemical  and  electronic  structure. 

4.  RESULTS  AND  DISCUSSION 

4.1.  Film  Growth  and  Characterization 

All  amorphous  hydrogenated  boron  carbide  (a-BvC:Hv)  films  described  in  this  report  were 
grown  by  plasma-enhanced  chemical  vapor  deposition  (PECVD)  with  an  ortho- carborane  (o- 
C2B10H12)  precursor  and  argon  process  gas  using  a  custom-built  capacitively  coupled  system. 
Details  of  the  thin-film  deposition  process  are  described  in  Nordell  et  al.  JAP  118  (2015) 
035703.23  Briefly,  the  system  consists  of  a  12  cm  diameter  upper  rotating  anode  and  substrate 
holder/heater  and  a  25  cm  diameter  lower  cathode  and  precursor  delivery  system  (i.e.,  showerhead) 
housed  within  an  80  L  chamber.  RF  power  at  a  standard  frequency  of  13.56  MHz  is  applied  across 
the  plates,  at  a  separation  of  4.4  cm.  To  deliver  the  solid-state  precursor  to  the  showerhead,  heated 
argon  gas  is  flowed  through  a  heated  ‘bubbler’  system,  which  consists  of  a  glass  nipple  containing 
ortho-c arborane  loaded  on  top  of  a  stack  of  glass  beads.  The  bubbler  is  heated  to  75  °C,  and  the 
gas  lines  and  showerhead  to  90  °C.  The  ortho-c  arborane  is  sourced  from  Katchem,  and  resublimed 
in  vacuo  prior  to  use.  The  argon  process  gas  is  sourced  from  Airgas  at  a  purity  of  <10  ppb  O2  and 
<20  ppb  H2O.  Restek  O2  and  H2O  filters  are  installed  in  line  to  reduce  the  levels  of  these 
contaminants  to  the  ppt  range.  Typically,  batches  of  films  are  grown  on  different  substrates  so  that 
these  may  be  characterized  by  multiple  techniques  including  ellipsometry  (thickness,  optical 
constants),  optical  transmission  spectroscopy  (absorption  coefficient,  band  gap,  Urbach  energy, 
Tauc  parameter),  photoemission  spectroscopy  (atomic  composition,  electronic  structure),  fourier 
transform  infrared  spectroscopy  (chemical  bonding),  nuclear  reaction  analysis  (atomic 
composition,  including  at.%  H),  nanoindentation  (hardness,  Young’s  modulus),  and  current- 
voltage  and  capacitance-voltage  measurements. 
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We  have  designed  a  number  of  film  growth  experiments  to  systematically  test  the  influence 
of  process  parameters  on  thin-film  properties.  These  are  summarized  in  Table  1.  Results  from  the 
B,  C,  D,  and  N  series  growths  are  described  in  a  series  of  publications: 

•  B.  J.  Nordell,  S.  Karki,  T.  D.  Nguyen,  P.  Rulis,  A.  N.  Caruso,  S.  S.  Purohit,  H.  Li,  S.  W.  King, 
D.  Dutta,  D.  Gidley,  W.  A.  Lanford,  and  M.  M.  Paquette,  ‘The  Influence  of  Hydrogen  on  the 
Chemical,  Mechanical,  Optical/Electronic,  and  Electrical  Transport  Properties  of  Amorphous 
Hydrogenated  Boron  Carbide,’  Journal  of  Applied  Physics  118  (2015)  035703. 

rhttp://dx.doi.org/10. 1063/1 .49270371 

•  B.  J.  Nordell,  C.  L.  Keck,  T.  D.  Nguyen,  A.  N.  Caruso,  S.  S.  Purohit,  W.  A.  Lanford,  D.  Dutta, 
D.  Gidley,  P.  Henry,  S.  W.  King,  and  M.  M.  Paquette,  ‘Tuning  the  Properties  of  a  Complex 
Disordered  Material:  Full  Factorial  Investigation  of  PECVD-Grown  Amorphous 
Hydrogenated  Boron  Carbide,’  Submitted  to  Materials  Chemistry  and  Physics,  August  19 
2015,  MATCHEMPHY S-D-15-01 944. 

•  M.  M.  Paquette,  J.  W.  Otto,  W.  A.  Lanford,  B.  J.  Nordell,  S.  Karki,  S.  W.  King,  and  A.  N. 
Caruso,  ‘Thermal  Treatment  and  Depth  Profile  Studies  of  Thin-Film  Carborane-Based 
Amorphous  Hydrogenated  Boron  Carbide  Growth  by  Plasma -Enhanced  Chemical  Vapor 
Deposition,’  (2015),  To  be  submitted. 


Table  1.  A  selection  of  PECVD  a-BrC:Hv  film  growth  experiments. 


Series 

No 

films 

Description 

Growth 
temp  (C) 

RF  Power 
(W) 

Pressure 

(Torr) 

Total  flow 
rate  (seem) 

Partial 
flow  rate 

Anneal 

conditions 

B 

5 

Temp 

dependence 

300-500 

30 

200 

50 

0.2 

— 

C 

7 

Anneal 

dependence 

350 

30 

200 

50 

0.2 

Variable 

D 

24 

Power/temp 

25-450 

5-30 

200 

50 

0.2 

- 

N 

37 

25  factorial 

50,  400 

10,  40 

0.2,4 

50, 100 

0.1,  0.5 

- 

Q 

45 

One  at  a  time 

75-400 

10-40 

0.4-8 

10-200 

0.1-2* 

Variable 

p 

10 

Anneal 

400 

40 

0.4 

100 

0.5 

Variable 

dependence 
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Low  temp,  power: 

•  Higher  at.%  H  (>40%) 

•  Lower  density  (<1  g/cm3) 

•  Higher  pore  volume  (0.7  nm) 

•  Lower  hardness  (<2  GPa)  & 


High  temp,  power: 

•  Lower  at.  %H  (<1 5%) 

•  Higher  density  (>2  g/cm3) 

•  No  measurable  pores 


•  Higher  hardness  (>20  GPa)  & 
Young’s  modulus  (>300  GPa) 
•  Lower  £g  (<2.0  eV) 


•  Higher  k(>6) 


Figure  1.  Atomic  concentration  hydrogen  (at.%  H)  in  a-BvC:Hv  films  as  a  function  of  varying 
temperature  and  R.F.  power  during  PECVD  growth  demonstrating  how  a-BvC:Hv  can  range  from 
low-density  hydrogen-rich  to  high-density  hydrogen-deficient  films.  [Figure  adapted  from 

NordeUet  al.  .TAP  118  (2015)  03570323!. 

The  first  publication,  which  covers  a  series  of  a-B.vC:Hv  films  deposited  with  varying  growth 
temperature  and  power  (i.e.,  the  “D”  series),  emphasizes  the  influence  of  hydrogen  concentration 
on  thin-film  properties.  We  found  that  we  were  able  to  produce  films  with  an  extremely  wide  range 
of  properties,  including  H  concentration  from  10-45%,  density  from  0.9-2. 3  g/cm3,  hardness  from 
1-25  GPa,  Young’s  modulus  from  10-340  GPa,  Tauc  optical  band  gap  from  1.7-3. 8  eV,  Tauc 
slope  parameter  (BU1)  from  -200-1600  cm_1/2-eV_1/2,  Urbach  energy  (£u)  from  0. 1-0.7  eV, 
dielectric  constant  from  3. 2-7. 6,  and  electrical  resistivity  from  1010— 1015  Q-cm  (Figure  1). 
Hydrogen  concentration  was  found  to  correlate  directly  with  thin-film  density,  and  both  were  used 
to  map  and  explain  the  other  material  properties  (Figure  2).  Hardness  and  Young’s  modulus  exhibit 
a  direct  power  law  relationship  with  density  above  -1.3  g/cm3  (or  below  -35%  H),  below  which 
they  plateau,  providing  evidence  for  a  rigidity  percolation  threshold.  An  increase  in  band  gap  and 
decrease  in  dielectric  constant  with  increasing  H  concentration  were  explained  by  a  decrease  in 
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network  connectivity  as  well  as  mass/electron  density.  An  increase  in  disorder,  as  measured  by  the 
parameters  Eu  and  BV2,  with  increasing  H  concentration  was  explained  by  the  release  of  strain  in 
the  network  and  associated  decrease  in  structural  disorder.  All  of  these  correlations  in  a-BAC:Hv 
were  found  to  be  very  similar  to  those  observed  in  amorphous  hydrogenated  silicon  (a-Si:H),  which 
suggests  parallels  between  the  influence  of  hydrogenation  on  their  material  properties  and  possible 
avenues  for  optimization. 


(i) 


(b) 


(f) 


(c) 


(g) 


(d) 


(h) 


0) 


(k) 


Figure  2.  Correlations  between  a  series  of  material  properties  and  atomic  concentration 
hydrogen  in  a-BAC:Hy  films:  (a)  density,  (b)  pore  diameter,  (c)  hardness  (H),  (d)  Young’s 
modulus  (E),  (e)  Tauc  optical  band  gap  (iiTauc,  n  =  2),  (f)  Eq4  optical  band  gap,  (g)  Urbach  energy 
(Fu),  (h)  Tauc  slope  parameter  (Bu2,  n  =  2),  (i)  high-frequency  dielectric  constant  (ei),  (j)  total 
(low-frequency)  dielectric  constant  (a;),  (k)  k  -  ei,  and  (1)  electrical  resistivity  (p).  [Figure  from 

Nordell  et  al.  .TAP  118  (2015)  03570323!. 
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The  second  publication  is  based  on  a  25  full  factorial  experiment.  This  is  an  experiment  in 
which  all  unique  combinations  of  each  of  five  process  parameters  at  two  levels,  “low”  and  “high,” 
are  tested,  which  allows  for  a  statistically  sound  evaluation  of  both  the  effect  of  each  parameter  as 
well  as  any  interactions  between  parameters,  while  minimizing  the  number  of  growths  required. 
This  study  expanded  upon  the  previous  work  looking  at  growth  temperature  and  power,  to  include 
the  effects  of  other  major  PECVD  process  parameters:  pressure,  total  flow  rate,  and  partial 
precursor  flow  rate.  Ultimately,  the  films  produced  demonstrated  a  similar  range  of  properties; 
however,  from  the  results  we  were  able  to  extract  a  more  thorough  understanding  of  the  effects  of 
process  parameters  on  thin-film  properties. 

One  of  the  important  results  highlighted  in  this  publication  was  the  wide  range  of  growth  rates 
possible,  which  extended  from  1-112  (±5)  nm/min.  Pressure  and  partial  precursor  flow  rate  were 
found  to  have  the  greatest  effect  on  growth  rate,  with  low  pressure  and  high  partial  flow  rate 
conditions  yielding  the  highest  growth  rates;  this  can  be  explained  by  the  higher  delivery  rate  of 
reactive  species  to  the  substrate  surface  under  these  conditions.  Power  and  total  flow  rate  also 
demonstrated  smaller,  yet  still  statistically  significant,  effects.  Temperature  did  not  show  any 
effect,  which — combined  with  the  above  results — suggests  that  growth  rate  is  primarily  limited  by 
mass  transport  rather  than  reaction  rate.  These  results  are  important  in  the  context  of  neutron 
detection,  since  higher  growth  rates  are  desirable  toward  producing  relatively  thick  (50-100  pm 
vs  <1  pm)  a-BYC:Hv  films. 

A  second  important  result  concerns  the  influence  of  growth  conditions  on  thin-film  properties. 
Hydrogen  concentration  and  density  were  found  to  be  influenced  primarily  by  temperature  and 
power,  consistent  with  the  previous  study.  However,  the  optical  and  electronic  structure  properties 
were  found  to  be  influenced  not  only  by  temperature,  but  also  by  a  significant  power*pressure 
interaction  effect.  As  a  result,  the  direction  of  the  power  effect  depends  on  whether  one  is  operating 
at  high  or  low  pressure  and  vice  versa.  Thus,  two  opposing  growth  regimes  can  be  identified — 
high  power  and  high  pressure  or  low  power  and  low  pressure — which  are  both  able  to  produce 
films  with  similar  optical  and  electronic  properties. 
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Figure  3.  A  series  of  properties  in  a-BAC:H,  films  vs  atomic  concentration  hydrogen  and  B/C 
ratio:  (a)  growth  rate,  (b)  high-frequency  dielectric  constant,  ei,  (c)  low-frequency  (total) 
dielectric  constant,  k,  (d)  band  gap,  £04,  (e)  Urbach  energy,  Ev,  (f)  Tauc  parameter,  Bul,  (g) 
electrical  resistivity  (log  p),  (h)  field-dependent  mobility  (log  uy),  and  (i)  charge  carrier 
concentration  (log  n).  [Figure  from  Nordell  et  al.  (2015)  submitted,  MATCHEMPHYS-D-15- 

019441. 


A  third  result  deals  with  a  deeper  look  at  carbon  concentration,  which  was  more  variable  in 
this  second  study.  The  metric  reported,  B/C  ratio,  was  influenced  by  temperature,  power,  pressure, 
and  the  power*pressure  interaction.  Interestingly,  it  was  found  that  B/C  ratio  correlated  strongly 
to  growth  rate,  but  was  completely  independent  from  hydrogen  concentration.  Through  an  analysis 
of  correlations  between  material  properties  (optical,  electronic,  electrical)  and  chemical 
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composition,  it  was  found  that  the  majority  correlated  strongly  to  atomic  concentration  hydrogen 
and  weakly-to-moderately  to  B/C  ratio  (Figure  3).  Thus,  although  H  concentration  (i.e.,  network 
connectivity/density)  appears  to  be  the  strongest  determiner  of  electronic  and  optical  properties, 
the  C  concentration  still  has  some  effect.  In  fact,  all  of  the  optical  and  electronic  structure 
parameters  investigated  showed  the  best  correlation  to  a  combined  parameter,  the  product  of  at.% 
H  and  B/C  ratio.  Unexpectedly,  the  charge  transport  properties  did  not  show  a  clear  correlation 
with  electronic  structure,  which  suggests  that  the  mechanism  cannot  be  explained  by  a  simple 
model  such  as  multiple  trapping.  This  will  be  discussed  further  in  Section  4.4. 


Etch  Depth  /  nm  Etch  Depth  /  nm 

Figure  4.  (a)  X-ray  photoemission  spectroscopy  depth  profile  data  for  a  low-density  film 
(B4.3COo.9:H4.8,  1.1  g/cm3).  (b)  B/C  ratio  as  a  function  of  etch  depth  obtained  from  data  in  (a). 


The  third  publication  describes  the  effect  of  annealing  time  on  B/C  ratio,  at.%  H  and  density 
in  a-BAC:Hv  films,  where  it  is  found  that  these  properties  continue  to  change  over  long  annealing 
times  (multiple  hours)  when  annealed  at  350  °C.  As  such,  the  at.%  H  and  at.%  C  decrease,  while 
the  density  increases  Additionally,  this  publication  describes  the  results  of  an  XPS  depth  profile 
study  on  a-B.vC:H_v  films  (Figure  4),  which  leads  to  two  main  conclusions:  (1)  that  the  oxide  layer 
is  constrained  to  the  top  few  nm’s  of  the  film,  even  after  the  prolonged  atmospheric  exposure  of  a 
relatively  low-density  film,  and  (2)  that  an  extended  carbon  gradient  exists,  where  the  surface  is 
carbon  rich  with  stoichiometry  of  ~a-B3C:Hv,  and  a  stable  “bulk”  composition  is  not  achieved  until 
>100  nm  into  the  film.  The  consequence  of  (2)  is  that  the  surface  and  bulk  properties,  which  are 
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known  to  correlate  to  chemical  composition,  may  differ,  and  that  surface-sensitive  characterization 
techniques  (e.g.,  XPS),  may  not  accurately  reflect  the  character  of  the  bulk. 

4.2.  Metal-to-a-BC:H  Interface  Studies 

4.2.1.  Cr,  Ti,  At,  and  Cu 

The  electronic  and  chemical  structure  of  the  metal-to-a-BxC:Hy  interface  was  studied  using 
valence  band  and  core  level  X-ray  and  ultraviolet  photoemission  spectroscopies  (XPS  and  UPS) 
for  a  series  of  low-work-function  metals  (Cr,  Ti,  Al,  and  Cu).  These  were  chosen  as  they  had 
comparable  work  functions  to  a-B*C:Hv,  with  Cr  in  particular  being  a  metal  that  had  been  used 
previously  as  a  contact  in  studies  of  a-BYC:Hv  neutron  detector  devices.  Films  were  grown  as 
described  in  Nordell  et  al.  .TAP  118  (2015)  03570323  with  a  growth  temperature  of  350  °C,  power 
of  30  W,  pressure  of  200  mTorr,  total  flow  rate  of  36  seem,  and  partial  flow  rate  of  0.5.  These 
growths  were  not  deposited  under  ultra-strict  oxygen  mitigation  conditions,  and  used  research 
grade  Ar  gas  without  additional  filters.  The  thermally  treated  sample  was  heated  at  475  °C  for  60 
min.  XPS  and  UPS  measurements  were  first  done  on  the  a-BvC:Hv  films  to  determine  atomic 
composition,  work  function,  and  valence  band  character.  The  films  measured  were  found  to  have 
a  composition  of  ~a-B3Co.5:Hv  (-10%  oxygen),  work  function  of  4.7  eV,  and  valence  band  maxima 
of  0.8  eV.  A  band  gap  value  of  2.4  eV  was  obtained  via  a  Tauc  analysis  of  optical  absorption 
spectra.  Two  main  studies  were  performed  to  investigate  the  extent  of  interfacial  diffusion,  new 
chemical  species  formed,  and  interfacial  electronic  structure  (Figure  5). 

(i)  A  metal  overlayer  was  evaporated  onto  a  sputter-cleaned  a-BAC:Hv  film  to  a  given  thickness, 
and  core  level  and  valence  band  photoemission  spectra,  sensitive  to  the  top  several  nm’s  of  the 
sample,  were  subsequently  acquired.  This  process  was  repeated  for  metal  thicknesses  of  0.5, 
1,2,5,  10,  150,  and  -250  A. 

o 

(ii)  A  5  A  thick  metal  overlayer  was  evaporated  onto  a  sputter-cleaned  a-BTC:Hv  film,  and  core 
level  and  valence  band  photoemission  spectra  were  acquired  first  for  the  sample  as-is,  and  then 
after  6,  12,  and  18  min  of  Ar+  etching. 
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Figure  5.  (a)  Schematic  of  photoemission  experiment  conducted  to  study  the  metal-to-a-B.vC:Hy 
interface,  and  (b)  representation  of  band  bending  analysis  by  monitoring  core  level  shift. 


Full  details  for  these  measurements  and  analyses  are  described  in  Driver  et  al.  JPCM  24 
(2012)  445001  24  The  major  conclusion  was  that  Cr,  Ti,  and  Al  all  formed  interfacial  oxides.  These 
are  anticipated  to  perturb  the  injection  barrier  and  the  local  electric  field  at  the  interface.  Cu, 
however,  forms  a  clean  metallic  interface.  Examples  of  metal  core  level  spectral  and  valence  band 
spectra  for  Cr  and  Cu  are  given  in  Figure  6  and  Figure  7,  respectively.  Tentative  band  diagrams 
for  each  of  the  interfaces  were  proposed  (Figure  8).  Band  bending  was  monitored  by  following  the 
shift  in  B  Is  core  level  energy  as  a  function  of  metal  overlayer  coverage,  and  the  shifts  were 
generally  consistent  with  those  predicted  for  the  metal  or  metal  oxide  interfaces. 
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Figure  6.  Cr  2p  core  level  photoemission  spectra  [(a)  as-deposited  and  (b)  thermally  treated]  and 
valence  band  photoemission  spectra  [(c)  as-deposited  and  (d)  thermally  treated]  for  a-BiCOo.siHv 
films  with  Cr  overlayers  as  a  function  of  overlayer  thickness.  [Figures  from  Driver  et  al.  JPCM 

24  (2012)  44500124]. 
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Figure  7.  Cu  2p  core  level  photoemission  spectra  [(a)  as-deposited  and  (b)  thermally  treated]  and 
valence  band  photoemission  spectra  [(c)  as-deposited  and  (d)  thermally  treated]  for  a-BsCOo  siHy 
films  with  Cu  overlayers  as  a  function  of  overlayer  thickness.  [Figures  from  Driver  et  al.  JPCM 

24  (2012)  44500124]. 
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Figure  8.  Band  diagrams  for  the  metal-to-a-B3COo.5:Hv  interface  for  Cr  (a),  Ti  (b),  Al  (c),  and 
Cu  (d)  taking  into  account  interfacial  oxide  layers  observed  by  XPS.  [Figure  from  Driver  et  al. 

JPCM  24  (2012)  44500124!. 

4.2.2.  Ni  and  Pd 

Based  on  the  results  from  the  first  series  of  metals,  additional  metals  were  chosen  that  had 
higher  work  functions  such  that  a  lower  barrier  to  hole  injection  and  higher  barrier  to  electron 
injection  could  be  achieved.  Specifically,  we  elected  to  test  Ni  (work  function  of  5.15  eV)  and  Pd 
(work  function  of  5.12  eV).25  Metal  overlayer  experiments  were  performed  in  a  similar  fashion  as 
for  the  previous  series  of  metals.  The  overlayers  of  Ni  were  applied  to  a  film  of  composition  a- 
B4.iCOo.3No.2:Hy  (i.e.,  6  at.%  O  and  3  at.%  N),  which  had  been  Ar-ion-etched  for  2  h  to  remove 
-20  nm  of  material.  The  Ni  2p  core  level  spectra  were  deconvoluted  into  six  peaks  to  model  its 
complex  satellite  structure  (Figure  9),  as  per  Grosvenor  et  al.26  The  spectrum  at  10  A  coverage 
appears  identical  to  that  of  metallic  Ni,  both  in  terms  of  energy  and  peak  envelope,  with  the  main 
peak  falling  at  852.6  eV,  identical  to  the  literature  value.26-27  Reported  energies  and  peak  envelopes 
for  nickel  oxides  are  quite  distinct  from  those  of  metallic  Ni.28-29  At  lower  coverages,  the  Ni  2p 
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o 

core  level  peak  is  shifted  to  slightly  higher  energies  (e.g.,  853.4  eV  at  2  A  coverage).  However,  no 
evidence  for  the  presence  of  NiO,  with  an  expected  peak  maximum  at  854.0  eV  and  a  prominent 
shoulder  at  -856  eV,  is  observed.  In  addition,  no  evidence  for  NiO  is  observed  from  the  O  Is 
spectrum  in  the  form  of  a  diagnostic  peak  at  529.6  eV;  the  spectrum  looks  identical  to  that  typically 
observed  for  a-BAC:Hv  films  with  a  small  amount  of  oxygen  contamination. 


o 

Figure  9.  Ni  2p  core  level  photoemission  spectra  for  an  a-B4.iCOo.3No.2:Hv  film  with  2  A  (a)  and 

o 

10  A  (b)  Ni  overlayers. 


Binding  Energy  /  eV 


o 

Figure  10.  Valence  band  photoemission  spectrum  for  an  a-B4.iCOo.3No.2:Hv  film  with  0  A 

o  o 

(bottom  trace,  black),  5  A  (middle  trace,  red),  and  10  A  (top  trace,  blue)  of  Ni. 
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Figure  11.  (a)  Band  diagram  for  a  clean  Ni-to-a-B4.iCOo.3No.2:Hv  interface,  (b)  B  Is  core  level 

o 

photoemission  spectrum  displaying  shift  in  peak  maximum  from  188.5  eV  at  0  A  coverage  to 

188.2  eV  at  10  A  coverage. 


The  valence  band  spectrum  (Figure  10)  is  also  diagnostic  of  metallic  Ni  being  formed 
immediately  at  the  a-B4.iCOo.3No.2:Hrto-nickel  interface.30  A  tentative  band  diagram  for  a  clean 
a-B4.iCOo.3No.2:Hy-to-Ni  interface  is  given  in  Figure  11(a).  Because  the  work  function  of  Ni  (5.15 
eV)  is  slightly  higher  than  that  of  amorphous  hydrogenated  boron  carbide  (4.7  eV),  an 
accumulation-type  contact  is  expected  (desirable  for  hole  transport  from  the  semiconductor  to  the 
metal)  with  a  low  Schottky  barrier.  A  small  shift  in  the  energy  of  the  B  Is  core  level  from  188.5 

o  o 

eV  at  0  A  coverage  to  188.2  eV  at  10  A  coverage  [(Figure  11(b)]  is  consistent  with  the  expected 
direction  of  band  bending,  and  therefore  supports  this  model. 

The  relative  atomic  concentrations  of  B,  C,  O,  N,  and  Pd  in  the  sampled  layer  for  an  a- 
B4.iCOo.3No.2:Hv  film  with  palladium  overlayers  were  determined.  Because  the  Pd  3p3/2  and  O  Is 
core  level  peaks  overlap,  the  quantification  of  oxygen  requires  the  deconvolution  of  the  two; 
palladium  may  be  quantified  via  the  Pd  3d  level.  The  results  are  shown  in  Figure  12.  As  Pd 
overlayers  are  applied,  the  concentration  of  Pd  increases  while  that  of  B  decreases,  as  expected, 
but  importantly,  the  concentration  of  O  increases  relative  to  that  in  the  bulk.  This  behavior  is 
similar  to  that  observed  in  the  cases  of  Cr,  Ti,  and  A1  overlayers  on  a-B.vC:Hv,  which  formed 
interfacial  metallic  oxides,  suggesting  that  a  similar  phenomenon  is  occurring  for  Pd.  Indeed,  in 

o 

the  Pd  3d  core  level  spectrum  [Figure  13(a)  and  (c)],  at  0.5  A  coverage,  palladium  monoxide  (PdO, 
Pd2+)  spectral  contributions  make  up  -85%  of  the  total  peak  area,  with  another  -10%  assigned  to 
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PdC>2  (Pd4+),  and  <5%  to  metallic  Pd0.31  At  higher  coverages,  an  increased  contribution  from 

o 

metallic  Pd  is  observed  in  the  Pd  3d  spectrum,  reaching  -65%  by  20  A  coverage  [the  remaining 
-35%  originating  from  PdO  (-30%)  and  PdCF  (<5%)].  An  accurate  analysis  of  contributions  to 
the  Ols  core  level  spectrum  is  impeded  by  the  existence  of  at  least  five  different  components  at 
nearly  the  same  energy  (O  Is  B-O,  O  Is  C-O,  O  Is  PdO,  Pd  3p3/2  Pd°,  and  Pd  3p3/2  Pd2+). 
Interestingly,  in  the  B  Is  spectrum  [Figure  13(b)  and  (d)],  as  overlayers  are  applied,  a  distinct  peak 
at  -193  eV  appears.  While  it  is  common  to  observe  a  peak  at  190-191  eV  in  the  B  Is  spectrum  of 
a-B.vC:Hv  films  with  oxygen  contaminants,  which  we  have  assigned  to  boron  suboxides  (B1+/B2+), 
the  peak  at  193  eV  has  not  previously  been  observed  for  any  of  the  a-BCO:H  films  we  have  grown, 
nor  at  the  a-BCO:H-to-metal  interface  for  any  of  the  metals  studied  (including  Cr,  Ti,  Al,  Cu,  and 
Ni)  and  can  be  assigned  to  boron  oxides  in  the  3+  oxidation  state,  namely  B2O3  or  H3BO3. 32,33  It 
is  possible  that  Pd  uniquely  catalyzes  a  boron  oxidation  reaction. 


Overlayer  Thickness  /  Angstrom 

Figure  12.  Relative  atomic  concentrations  of  B,  C,  O,  N,  and  Pd  (at.%)  in  the  sampled  layer  (i.e., 
the  top  several  nm’s  of  the  sample)  determined  from  normalized  XPS  core  level  spectral 
intensities  (for  the  B  Is,  C  Is,  O  Is,  N  Is,  and  Pd  3d  core  levels)  for  an  a-B^iCOcuNo^FFy  film 
with  Pd  overlayers  as  a  function  of  overlayer  thickness. 
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Figure  13.  Pd  3d  and  B  Is  core  level  photoemission  spectra  for  an  a-B4.iCOo.3No.2:Hv  film  with 
Pd  overlayers  at  0.5  A  coverage  [(a)  and  (b),  respectively]  and  at  20  A  coverage  [(c)  and  (d), 

respectively] . 

The  formation  of  PdO  at  the  immediate  a-B4.iCOo.3No.2:Hv-to-Pd  interface  is  corroborated  by 
the  valence  band  spectrum  (Figure  14),  wherein  the  maximum  for  the  Pd  3d  peak  shifts  from  ~2.6 

o  o 

eV  at  0.5  A  coverage,  diagnostic  of  the  formation  of  PdO  at  low  coverages,  to  ~1.6  eV  at  20  A 
coverage,  diagnostic  of  metallic  Pd  at  thicker  coverages.34  35  A  tentative  band  diagram  for  an  a- 
BCO:H-to-Pd  interface  with  an  interfacial  PdO  layer  is  given  in  Figure  15.  Little  information  is 
available  in  the  literature  regarding  the  electronic  structure  of  PdO.  A  work  function  of  5.9  eV  and 
band  gap  of  -1  eV  were  used  to  construct  the  diagram,34,36  although  it  must  be  noted  that  these 
values  are  not  rigorously  established.  Assuming  they  are  adequate  estimates,  it  may  be  concluded 
that  the  PdO  layer  will  lead  to  significant  band  bending  and  space-charge  region  formation  in  the 
a-B.vC:Hv  layer.  It  is  possible,  given  that  the  PdO  layer  is  relatively  thin  (1-2  nm)  and  that  PdO  is 
a  narrow-band-gap  semiconductor,  that  its  presence  will  not  significantly  impede  charge  transport 
at  the  junction.  As  regards  an  experimental  determination  of  band  bending  through  photoemission 
spectroscopy,  no  shift  was  observed  for  the  Bis  core  level,  which  suggests  that  band  bending  may 
be  much  less  than  implied  in  Figure  15. 
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Figure  14.  Valence  band  photoemission  spectrum  for  an  a-B4.iCOo.3No.2:Hy  film  with  Pd 
overlayers  as  a  function  of  overlayer  thickness  (0,  0.5,  1,  2,  5,  10,  20  A). 


Figure  15.  Band  diagram  for  a  Pd-to-a-B4.iCOo.3No.2:Hy  interface  with  an  interfacial  PdO  layer. 


In  summary,  Ni  appears  to  form  a  clean  interface  to  a-EUiCOo.yNo^Hy,  with  no  evidence  for 
the  formation  of  new  chemical  species.  On  the  basis  of  band  alignment  predictions,  and  consistent 
with  experimentally  observed  band  bending,  Ni  would  be  expected  to  make  a  good  ohmic  contact 
to  a-BvCrHy,  similarly  to  Cu.  Palladium,  on  the  other  hand,  forms  a  thin  interfacial  oxide  layer, 
made  up  predominantly  of  PdO.  A  band  diagram  for  the  oxidized  Pd-to-a-BvC:H}  interface  is 
proposed,  though  it  is  uncertain  whether  this  layer  will  prove  advantageous  or  disadvantageous  for 
charge  transport  at  the  junction  due  to  the  scarcity  of  electronic  structure  information  for  PdO. 
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4.2.3.  Metal-Insulator-Metal  Heterostructure  Characterization 

Metal-insulator-metal  (MIM)  heterostructures  were  fabricated  to  determine  whether  a  given 
metal  formed  an  ohmic  junction  to  a-BYC:Hv.  This  may  be  evaluated  through  an  analysis  of  the  I- 
V  behavior  of  the  heterostructure,  which,  if  ohmic,  will  be  linear  (i.e.,  obey  Ohm’s  law,  I  =  V/R). 
MIM  heterostructures  were  fabricated  successfully  using  various  combinations  of  Cu,  Al,  Ti,  and 
Pd  contacts.  The  fabrication  of  Ni/a-BAC:Hv/Ni  heterostructures  was  unsuccessful  due  to  the  a- 
B.yC:Hv  films  consistently  delaminating  from  the  base  evaporated  Ni  layer.  The  fabrication  of  Cu/a- 
BAC:HV/Ti  heterostructures  was  also  unsuccessful  due  to  the  delamination  of  Ti  from  the  a-BAC:Hv. 
Additional  work  on  adhesion  layers  may  be  needed  to  further  develop  these  metals  as  contacts. 

Current-voltage  measurements  were  completed  in  ‘loops’  in  order  to  investigate  a  possible 
hysteresis  effect.  The  current  was  thus  measured  from  0  V  — >  x  V  — >  —x  V  — >  0  V,  where  x  was 
varied  from  low  to  high  voltage  (0.5,  1,  2,  5,  10,  15,  20  V)  until  erratic  behavior  was  observed.  For 
each  voltage  range,  three  full  cycles  were  measured  to  verify  electrical  stability. 

I-V  curves  for  a  series  of  devices  are  shown  in  Figure  16:  Cu/a-B.yC:Hy/Cu,  fabricated  entirely 
in  situ  (a);  Cu/a-BAC:Hv/Cu,  with  the  bottom  contact  exposed  to  atmosphere  (b),  Cu/a-BxC:Hy/Pd, 
with  the  bottom  contact  exposed  to  atmosphere  (c),  and  Al/a-B.yC:Hy/Al,  again  with  the  bottom 
contact  exposed  to  atmosphere  (d).  It  can  be  seen  that  the  devices  containing  different  metal 
contacts,  fabricated  both  in  situ  and  ex  situ,  display  similar  curve  shapes  and  current  magnitudes 
(within  the  same  order  of  magnitude);  the  differences  are  subtle.  All  of  the  curves  display 
hysteresis,  and  all  appear  to  display  linear  behavior  (/  cc  V)  at  low  voltages/fields,  and  curvature  (/ 
oc  V">1)  at  higher  voltages/fields  (see  Section  4.3.1  for  a  more  detailed  discussion  of  /-V  behavior). 
In  comparing  Cu/a-BAC:Hv/Cu  devices  fabricated  in  situ  and  ex  situ  [(a)  and  (b)],  the  devices 
exhibit  very  similar  behavior;  and  in  fact,  the  device  prepared  partially  ex  situ  appears  to  display 
‘cleaner’  behavior  (i.e.,  the  current  overlaps  at  each  voltage  range).  The  Al/a-BAC:H,/Al  device 
(d),  expected  to  possess  an  insulating  AI2O3  interfacial  layer  due  to  ex  situ  preparation,  also  shows 
very  similar  behavior  (except  that  it  breaks  down  above  10  V),  which  suggests  that  charge  can 
easily  tunnel  through  the  oxide  barrier.  We  have  additionally  tested  MIS  heterostructures,  with  Si 
as  the  semiconductor,  and  in  some  cases  a  temporary  mercury  contact  (using  a  Hg  probe)  as  the 
metallic  electrode.  A  comparison  of  I-V  traces  is  shown  in  Figure  17,  where  Cu/a-BAC:Hy/Cu,  Si/a- 
B.vC:Hv/Hg  (with  both  degenerate  and  1-15  Q-cm  Si),  and  Si/a-BAC:Hv/Cu  heterostructures  all 
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exhibit  relatively  similar  I-V  behavior.  Ohmic  behavior  has  also  been  observed  in  Cu/a- 
B.vC:Hy/Ti/Cu,  Cu/a-BAC:Hv/oxide/Pd  (a-BrC:Hy  exposed  to  atmosphere  before  deposition  of  top 
contact),  Pd/a-B.vC:Hv/Pd,  and  ITO  (indium  tin  oxide)/a-BAC:Hv/ITO  heterostructures. 


Figure  16.  I-V  curves  for  four  different  a-BAC:Hrbased  metal-insulator-metal  devices  showing 
the  effects  of  different  metals  and  in  situ  vs  ex  situ  processing. 
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Figure  17. 1-V  traces  for  four  different  a-BvC:Hv  MIM  heterostructures. 


Overall,  we  conclude  that  the  contact  type,  which  includes  the  possibility  of  a  metal  or  a 
semiconductor  (e.g.,  Si)  with  or  without  the  presence  of  an  interfacial  oxide  layer,  does  not  have  a 
significant  effect  on  the  general  1-V  behavior  of  the  devices  at  low/moderate  fields.  This  does  not 
necessarily  mean  that  the  contact  does  not  have  any  effect  on  the  device  performance.  At  very  high 
fields,  the  injection  mechanism  may  change  from  ohmic  to  schottky,  where  the  electric  field  leads 
to  an  increase  in  barrier  while  necessitating  greater  current  injection  to  sustain  ohmic-based 
transport.  The  reason  for  this  is  two-fold:  (1)  the  extremely  high  resistivity  of  a-BAC:Hv  implies 
that  very  little  charge  injection  is  needed  to  supply  sufficient  current  to  the  device.  Thus,  even  in 
the  presence  of  a  barrier,  ohmic  injection  is  sustained  with  very  little  charge  injection;  and  (2)  the 
amorphous  nature  of  a-BAC:Hv,  which  leads  to  the  presence  of  band  tails  consisting  of  localized 
tail  states  presents  the  possibility  of  charge  injection  via  the  more  accessible  mid-gap  states  rather 
than  via  the  conduction/valence  band  alone. 

4.2.4.  Summary 

The  tendency  to  form  metallic  oxides  at  the  metal-to-a-BAC:Hv  interface  suggests  that  oxygen 
levels  in  a-BAC:Hv  films  will  play  an  important  role  in  device  performance,  not  only  in  affecting 
the  bulk  a-BAC:Hv  properties,  but  also  in  modifying  heterostructure  interfaces.  The  incorporation 
of  oxygen  into  a-BAC:Hv  films  during  growth  is  ubiquitous,  consistent  with  the  strong  oxygen- 
gettering  tendencies  of  boron  carbide.  We  hypothesize  that  bulk  oxygen  will  remain  an  issue  even 
with  lower  oxygen  concentrations  (3-5%),  which  are  usually  still  present  in  films  with  the  lowest 
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oxygen  content.  If  a  non-oxygen-gettering  metal  contact  is  required,  Cu  or  Ni  are  viable  options. 
Although  interfacial  oxides  will  very  likely  have  some  sort  of  effect  on  heterostructure 
performance,  in  terms  of  gross  I—V  behavior,  these  appear  to  have  little  influence  on  the  low-field 
ohmic  injection  mechanism  observed  in  the  case  of  the  majority  of  the  MIM  heterostructures 
measured.  The  ease  of  forming  ohmic  contacts  for  a-BAC:Hv  has  been  rationalized  in  terms  of  its 
high  resistivity  and  the  presence  of  localized  mid-gap  states  arising  from  its  amorphous  structure. 

4.3.  Charge  Transport  Measurements 


4.3.1.  Mobility  Measurements 

Measuring  mobility  in  low-mobility  materials  is  a  nontrivial  task.37  Not  only  is  the  concept  of 
a  precise  mobility  value  nebulous  due  to  the  fact  that  mobility  is  dependent  on  experimental 
conditions  (e.g.,  device  thickness,  electric  field)  and  ultimately  on  the  distribution  of  charge 
carriers  within  the  density  of  states  over  the  course  of  the  experiment  (and  therefore  the  portion  of 
the  density  of  states  in  which  the  mobility  is  being  probed),38  but  experiments  devised  to  measure 
mobility  all  necessitate  that  a  different  specific  set  of  requirements  be  met  in  order  to  ensure 
accuracy — oftentimes  these  requirements  cannot  easily  be  tested  and  can  only  be  assumed.39  A 
recent  paper  by  Blakesley  et  al.  demonstrating  the  wide  variability  in  mobility  values  calculated 
for  the  same  materials  serves  as  a  cautionary  note.37 


Of  the  methods  available  for  measuring  mobility,  drift-mobility-based  methods  are 
particularly  suitable  for  high-resistivity  (low-mobility)  solids,  as  they  ultimately  scale  with  transit 
time,  which  increases  with  decreasing  mobility  and  thus  becomes  easier  to  measure.  In  contrast, 
in  Hall  methods,  the  Hall  voltage  decreases  with  decreasing  mobility,  and  renders  measurements 
increasingly  challenging.  By  measuring  the  transit  time,  ftr,  of  charge  carriers  through  a  material, 
mobility  can  be  extracted  through  the  relationship  of  ftr  to  drift  velocity,  v,  which  is  in  turn  related 
to  mobility,  as  per  Equation  1,  where  d  represents  film  thickness  and  E ,  electric  field. 


d  d 
v  iiE 


[1] 


While  the  light-based  time-of-flight  (TOF)  method  may  be  the  most  commonly  applied  drift- 
mobility  method,  the  transient  (or  dark- injection)  space-charge-limited  current  (DI-SCLC)  method 
presents  several  advantages  in  our  case,  including  the  use  of  injecting  contacts,  thinner  films,  and 
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greater  compatibility  with  dispersive  transport.  A  short  list  of  additional  mobility  determination 
techniques  for  high-resistivity/low-mobility  materials  includes  charge  extraction  by  linearly 
increasing  voltage  (CELIV),  impedance/admittance  spectroscopy,  field  effect  transistor  methods, 
and  steady-state  space-charge-limited  current  (SS-SCLC)  methods.39  Blakesley  et  al.37  argue  that 
for  benchmarking  purposes,  the  steady-state  space-charge-limited  current  method  is  the  optimum 
choice,  but  this  may  simply  be  due  to  its  relative  simplicity  and  ease  of  standardization  (precision) 
rather  than  its  accuracy.  WeiB  et  al.40  argue  that  the  corollary  DI-SCLC  method  is  superior  in  terms 
of  determining  absolute  values,  as  it  measures  the  drift  mobility  more  directly  via  a  direct  measure 
of  transit  time  and  does  not  involve  the  ambiguity  related  to  trapping  inherent  in  SS-SCLC 
analyses.  Ultimately,  however,  there  are  so  many  unknowns  (e.g.,  dominant  charge  carrier,  density 
of  states,  trap  profile,  etc.),  particularly  when  dealing  with  a  complex  material  such  as  a-BvC:Hv 
that  has  been  minimally  studied,  that  the  use  of  a  series  of  methods  should  provide  the  most 
information  regarding  the  validity  and  advantages/disadvantages  of  the  different  techniques.  Both 
steady-state  and  transient  SCLC  techniques  have  been  applied  to  a-BvC:Hv  samples. 

The  complete  theory,  results,  and  analysis  of  this  work  is  described  in  detail  in  the  following 
thesis: 

•  Christopher  L.  Keck,  ‘Electrical  Carrier  Mobility  Measurements  of  Amorphous 
Hydrogenated  Boron  Carbide  Using  Space-Charge-Limited  Current  Techniques,’  M.Sc. 
Thesis,  University  of  Missouri-Kansas  City  (2015). 

rhttps://mospace.umsystem.edu/xmlui/handle/10355/465451 

4.3. 1.1. Steady-State  Space-Charge-Limited  Current 

In  an  insulating  solid,  we  expect  to  observe  space-charge-limited  current  effects.42  When  an 
electric  field  is  applied  to  an  insulator  or  semi-insulator  sandwiched  between  conducting 
electrodes,  ohmic  behavior  (/  x  V)  will  typically  be  observed  at  low  field.  In  this  regime,  the 
contact  is  supplying  sufficient  charge  to  replace  that  within  the  solid.  In  this  scenario,  a  contact 
may  be  referred  to  as  ‘injecting’  or  ‘ohmic,’  which  does  not  necessarily  imply  good  band  lineup 
or  accumulation  layer  formation  (typical  markers  of  injecting/ohmic  contacts),  but  rather  a  low 
contact  resistance  relative  to  bulk  resistance,  which  is  a  much  easier  requirement  to  meet  in  the 
case  of  a  high-resistivity  solid.  As  the  electric  field  is  increased,  the  I—V  behavior  begins  to  take 
on  a  superlinear  relationship  with  I  ccV1  [Figure  1 8(a)] .  At  this  point,  injected  charge  accumulates 
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within  the  insulator  adjacent  to  the  electrode  faster  than  it  can  be  transported  through  the  device, 
causing  a  buildup  of  space  charge  and  modification  of  I—V dependence,  which  now  follows  Child’s 
law: 


9  V2 

Jsclc-  t2] 

where  s  =  coCi,  with  eo  representing  the  absolute  permittivity  of  free  space  and  cr  the  relative 
dielectric  constant  of  the  material.  The  distinct  transition  between  ohmic  and  SCLC  transport 
regimes  can  be  understood  as  the  field  at  which  dielectric  relaxation  time,  td  (=  s/a),  which 
effectively  allows  for  neutralization  of  excess  injected  charge,  becomes  longer  than  the  transit 
time,  Ztr  (=  d/v  =  d/juE)  (h  in  Figure),  which  becomes  shorter  with  increasing  electric  field 
proportionally  to  mobility  [Figure  18(b)].41  The  presence  ofthe  mobility  term  in  Child’s  law  allows 
for  a  determination  of  this  property  from  a  fitting  of  the  I—V  (or  J—V  or  J—E)  curve  within  the  SCLC 
regime.  The  downfall  of  this  method  is  that  Child’s  law  assumes  a  perfect  solid  with  no  trapping 
as  well  as  field-independent  mobility,  which  are  both  not  the  case  in  the  majority  of  insulators, 
particularly  those  exhibiting  disorder.  Corrections  to  the  expression  to  account  for  trapping  and  for 
field-dependent  mobility  exist,  but  traditional  treatments  are  still  oversimplifications.43  Aware  of 
the  limitations  of  the  model,  we  have  opted  to  apply  the  Murgatroyd44  treatment  which  takes  into 
account  the  field-dependent  mobility: 

V-f  =  Eo  exp (yVE)  [3] 


and  approximates  SCLC  current  density,  /sclc,  as: 


Jsclc 


9  V2  (  fv\ 
-^-explOMy^j 


[4] 


where  juo  represents  the  limit  of  mobility  at  zero  electric  field  and  y  describes  its  field-dependence. 
While  the  field  dependence  cannot  account  for  all  trapping  profiles,  it  can  approximate  shallow 
trapping  effects.  We  have  applied  the  protocol  outlined  by  Blakesley  et  al.37  in  the  hopes  of,  if  not 
obtaining  100%  accurate  mobility  values,  at  least  obtaining  consistent  and  precise  mobility  values 
that  can  be  compared  to  each  other,  and  to  others  reported  in  the  literature  obtained  by  a  similar 
protocol.  It  is  reasonable  to  assume  that  the  mobility  thus  obtained  can  represent  a  lower  limit  on 
this  property.40 
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Figure  18.  (a)  Ohmic  and  SCLC  I—V regimes;  and  (b)  distinction  between  regimes  as  a  function 
of  transit  time  and  dielectric  relation  time  constant  (adapted  from  ref41). 


Steady-state  current-voltage  traces  were  obtained  on  two  types  of  heterostructures, 
Cu/aBAC:Hv/Cu  and  Si/a-BAC:Hv/Hg  (Hg  represents  a  temporary  mercury  contact),  which  both 
demonstrate  very  similar  I—V  curves.  The  advantage  of  Hg  probe  measurements  is  the  ease  of 
fabricating  simple  thin-film  Si/a-BAC:Hv  heterostructures  relative  to  MIM  sandwich  devices,  and 
therefore  the  possibility  of  high  throughput  for  property  correlations.  Voltage  is  typically  varied 
cyclically  from  0  V  — >  x  V  — >  —x  V  — >  0  V,  while  increasing  voltage  limits,  x,  in  increments  until 
erratic  behavior  is  observed,  and  repeating  each  cycle  three  times  to  allow  the  device  to  stabilize 
and  ensure  reproducibility.  As  observed  in  Figure  19,  the  I-V  curves  exhibit  hysteresis, 
hypothesized  to  be  caused  by  charge  trapping/detrapping  which  induces  a  displacement  current 
due  to  changing  capacitance  in  addition  to  the  expected  resistive  current.45-47  The  nonlinearity  of 
the  I-V  curves  is  consistent  with  SCLC  effects:15  a  log-log  plot  of  the  data  (Figure  20)  highlights 
two  relatively  linear  regions  with  I  cc  V”  with  n  =  1  (ohmic)  and  n  =  2  (SCLC). 
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Figure  19.  Example  IV  plot  of  an  a-BYC:Hv  MIM  test  device  displaying  hysteresis.  [Figure  from 

Keck  M.Sc.  Thesis  (2015)1. 


Figure  20.  Example  log-log  IV  plot  for  a  Hg/a-B.vC:Hv/Si  test  device  (N4)  demonstrating  two 
distinct  linear  regions.  A  linear  fit  to  the  low-field  region  yields  a  slope  of  0.95,  consistent  with 
ohmic  behavior.  A  linear  fit  to  the  higher-field  region  yields  a  slope  of  1.8,  consistent  with  SCLC 

behavior.  [Figure  from  Keck  M.Sc.  Thesis  (2015)1 . 
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1  2  3  4  5  6  7 

Voltage  (V) 

Figure  21.  Fit  of  the  SCLC  region  from  the  data  shown  in  Figure  20  to  Murgatroyd’s  expression. 

[Figure  from  Keck  M.Sc.  Thesis  (2015)1 . 


The  ohmic  region  of  the  curve  can  be  fit  to  obtain  resistivity  via  the  standard  expression,  where  p 
is  resistivity,  V  is  the  applied  voltage,  I  is  the  measured  current,  A  is  the  contact  area,  and  d  is  the 
sample  thickness: 


P 


VA 

Id 


[5] 


Fitting  the  SCLC  region  of  the  curve  to  the  Murgatroyd  expression  (via  the  Blakesley  protocol37) 
in  the  form: 


IscLc=A1V2exp{t1yTv)  [6] 


„  _9Ap0s 
1  8d3 


 0.89 ly 


[7] 

[8] 


yielded  parameters  A\  and  ti,  from  which  po  and  y  could  be  calculated.  Inserting  these  terms  into 
Equation  3  yields  the  field-dependent  mobilities. 
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Resistivity  (p ),  mobility  (p),  and  free  carrier  concentration  in)  (obtained  via  p  =  l/enp)  were 
measured  for  a  large  number  of  samples.  Because  p,  p,  and  n  are  related  via  p  =  1 1 e mi,  we  expect 
a  correlation  between  variables.  Indeed,  there  is  a  strong  linear  correlation  between  log(//F)  and 
log(p)  [Figure  22(a)]  as  well  as  a  definite  correlation  between  log(/n  j  and  login)  [Figure  22(b)] 
and  login)  and  log!/;)  [Figure  22(c)].  In  the  latter  two  cases,  however,  the  log(/vi  j-log(/z)  and 
log(n)-log(yo)  scatter  plot  appears  to  define  a  plane  rather  than  a  line.  Indeed,  when  all  three 
variables  are  plotted  against  one  another  in  3D  space,  this  plane  is  clearly  defined  [Figure  22(d)] 
by  four  extremes:  low  p  (10“14),  high  n  (1017)  and  high  p  (1016);  high  p  (10“5),  low  n  (1015),  and 
low  p  (108);  low  p  (10“13),  very  high  n  (1019),  and  moderate  p  (1013);  and  moderate  p  (10“9),  low  n 
(1014),  and  moderate  p  (1013).  Overall,  this  large  data  set  demonstrates  the  extreme  tunability  of  a- 
B.yCHv,  with  mobility  for  example  ranging  across  ten  orders  of  magnitude,  and  allows  us  to 
understand  the  combinations  of  charge  transport  properties  possible. 


(a)  (b) 


Figure  22.  Correlations  between  (a)  mobility  (//i  ,  at  0.1  MV/cm)  and  electrical  resistivity  (p),  (b) 
pi  and  charge  carrier  concentration  (n),  (c)  p  and  n,  and  (d)  all  three  variables  (all  plotted  as 

logarithmic  values). 
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4.3. 1.2. Dark-Injection  Transient  Space-Charge-Limited  Current 

In  the  DI-SCLC  method,  a  voltage  pulse  is  applied  to  a  sandwich-geometry  unipolar 
heterostructure  (i.e.,  the  material  is  either  majority  electron/hole  or  has  injecting  contacts  for  one 
type  of  carrier  and  blocking  contacts  for  the  other),  and  the  time-dependent  current  response 
recorded  (Figure  23).  The  voltage  pulse  generates  electron-hole  pairs  at  the  surface  of  the  device, 
effectively  creating  a  ‘sheet’  of  charge  carriers  that  induces  a  transient  current  response  in  the 
external  circuit  as  it  is  swept  across  the  device  under  an  applied  electric  field.  The  current  transient 
will  display  a  maximum  at  time  t  =  r,  which  is  related  to  the  transit  time  of  the  sheet  of  charge,  tv, 
by:48 

r  =  0.787ttr  [9] 


Figure  23.  Idealized  transient  current  response  from  a  DI-SCLC  experiment.  [Figure  adapted 

from  ref49] . 


One  of  the  challenges  to  implementing  the  DI-SCLC  method  for  a  high-resistivity  solid  such 
as  a-BvC:Hv  is  surmounting  RC  time  delay  issues.  To  convert  current  into  a  measurable  voltage, 
measuring  across  a  resistor  with  a  sufficiently  large  resistance  is  necessary.  The  result  is  a  voltage 
decay  that  masks  the  desired  signal  [Figure  24(a)  and  (b)].  While  Scott  et  al.50  originally  proposed 
to  rectify  this  with  a  bridge  circuit  [Figure  24(c)],  we  have  opted  to  implement  a  transimpedance 
amplifier  circuit  [Figure  24(d)],  more  recently  described  by  Szymanski  et  al.51  Here,  current  is 
directed  into  the  negative  input  of  an  op  amp.  A  resistor  is  connected  across  the  op  amp’s  negative 
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input  and  output  while  the  positive  input  is  grounded.  An  oscilloscope  measures  the  voltage  from 
the  op-amp’s  output.  By  placing  the  resistor  across  the  op-amp,  its  resistance  is  removed  from  the 
circuit.  This  reduces  the  series  resistance  of  the  circuit  and  results  in  a  smaller  RC  time  constant 
that  will  not  override  the  current  transient  peak. 


Figure  24.  (a)  simple  resistor  circuit;  (b)  effect  of  RC  voltage  delay  (dotted  line)  swamping  out 
desired  signal  (solid  line)  [figure  from  ref52];  (c)  bridge  circuit;  (d)  transimpedance  amplifier 

circuit. 


Transient  current  measurements  were  completed  on  a  subset  of  samples.  Not  all  samples 
yielded  usable  results,  typically  because  of  either  electrical  breakdown  or  insufficient  output 
current.  Representative  data  is  shown  in  Figure  25.  From  the  current  transients,  peak  time,  r,  was 
extracted  by  taking  the  derivative  of  a  polynomial  fit,  and  the  mobility  calculated  based  on  transit 
time  obtained  from  Equation  9.  Since  the  mobility  in  a-B<C:Hv  exhibits  a  field  dependence,  the 
zero-field  mobility  (juo)  and  gamma  factor  (y)  were  extracted  by  fitting  a  plot  of  /a  vs  Em  (Figure 
26).  Results  for  both  steady-state  and  transient  measurements  are  compiled  in  Table  2,  and 
illustrated  for  select  samples  in  Figure  27. 
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Time  (s) 


Figure  25.  DI-SCLC  current  transients  for  a  Hg/a-BAC:Hv/Si  heterostructure  (N16)  at  varying 
fields.  [Figure  from  Keck  M.Sc.  Thesis  (2015)1. 


Figure  26.  Mobility  vs  EU2  for  the  data  from  Figure  24.  [Figure  from  Keck  M.Sc.  Thesis 

(2015)1. 
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Table  2.  Comparison  of  mobility  metrics  obtained  from  transient  (dark-injection)  and  steady-state  mobility  measurements. 


DI-SCLC 

SS-SCLC 

7 

(cm/V) 1/2 

//o 

(cm2/V*s) 

flF 

(0.5  MV/cm) 
(cm2/V-s) 

fiF  range 
(cm2/[V-s]) 

E range 
(kV/cm) 

7 

(cm/V) 1/2 

//o 

(cmW*s) 

JUF 

(0.5  MV/cm) 
(cm2/V-s) 

fiF  range 
(cm2/V-s) 

E range 
(kV/cm) 

A 

-1.1  x  10“2 

2.4  x  IO"5 

3.2  x  10“8 

5.0  x  10“6-2.2  x 

io45 

23-56 

-3.2  x  10“3 

2.7  x  IO4 

2.8  x  104fl 

1.7  x  104— 1.3  x 

104 

23-78 

N16 

3.6  x  10~3 

5.3  x  104° 

3.2  x  10~9 

2.8  x  104-1.9x 

10“8 

244-976 

5.6  x  104 

7.5  x  IO41 

1.1  x  1040“ 

1.5  x  104O-1.8  x 

lO-io 

1460-2440 

N28b 

-8.6  x  10-4 

8.9  x  10~9 

5.8  x  IO4 

6.0  x  10~9-4.9  x 

10“9 

134-470 

-1.1  x  10“3 

8.3  x  IO4 

3.8  x  IO4 

6.2  x  104-4.1  x 

104 

67-393 

2N32 

1.9  x  10“2 

2.2  x  1043 

3.5  x  104° 

2.6  x  104°-6.8  x 

IO40 

134-175 

1.1  x  IO4 

1.3  x  1043 

2.8  x  1043 

1.9  x  1043-2.0x 

1043 

103-155 

N40 

-1.6x  10“2 

1.8  x  1045 

2.2  x  1041 

4.4  x  1041-3.1  x 

1041 

436-473 

-2.0  x  10“5 

7.4  x  10~9 

7.3  x  10“9 

7.3  x  10_9-7.3  x 

10~9 

303-909 

Q14 

-1.3  x  10“3 

6.0  x  1010 

3.1  x  104° 

2.7  x  104O-1.8x 

IO40 

339-847 

-1.6  x  104 

8.1  x  1044 

7.2  x  1044 

7.7  x  1044— 7.5  x 

1044 

76-254 

Q19 

-2.6  x  104 

3.7  x  1040 

3.1  x  104° 

3.1  x  104°-3.1  x 

IO40 

460-540 

-5.5  x  104 

9.6  x  1044 

6.5  x  1044 

7.4  x  1044-6.4  x 

1044 

230-550 

Q26 

-1.6  x  10“3 

1.5  x  10“9 

6.7  x  104° 

7.2  x  104°-6.0  x 

IO40 

216-333 

-2.2  x  IO-3 

4.5  x  1044 

9.5  x  1044 

2.5  x  10“14— 1 .3  x 

1044 

69-314 

Q28 

-4.6  x  10“3 

2.5  x  10“9 

4.2  x  104° 

7.3  x  1040-2.2  x 

IO40 

71-286 

3.8  x  IO4 

7.8  x  1042 

1.2  x  1041 

2.2  x  1041-2.9  x 

1041 

79-125 

"Electric  field  range  did  not  actually  extend  to  0.5  MV/cm 
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Figure  27.  Comparison  of  DI-SCLC  and  SS-SCLC  field-dependent  mobility  for  select  samples: 
(a)  A  (Cu/a-B.rC:Hy/Cu),  (b)  N16  (Hg/a-B.YC:Hy/Si),  (c)  N28b  (Hg/a-B,C:Hv/Si),  and  (d)  Q14 
(Hg/a-B.vC:Hy/Si).  [Figure  from  Keck  M.Sc.  Thesis  (2015)1. 


The  mobility  values  obtained  from  transient  measurements  were  found  to  be  higher  than  those 
obtained  from  steady-state  measurements,  by  between  one  and  four  orders  of  magnitude,  but  the 
trends  were  generally  the  same.  WeiB  et  al.40  obtained  a  similar  result,  and  ascribed  it  to  the  fact 
that,  while  transient  experiments  measure  transit  time  directly,  slow  trapping  can  have  a  significant 
effect  on  steady-state  measurements.  Obtaining  negative  values  for  y  is  also  an  unusual  result,  but 
may  be  associated  with  disorder  in  amorphous  solids.53 

Overall  we  conclude  that  steady-state  space-charge-limited  current  measurements  appear  to 
be  an  efficient  technique  for  high-throughput  mobility  measurements  on  a-BYC:Hv  thin  films,  as 
transient  measurements  were  found  to  be  fairly  demanding  to  implement  consistently  for  all  types 
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of  films  (e.g.,  different  thicknesses,  mobilities,  contact  types).  The  trends  in  mobility  values  were 
supported  by  both  steady-state  and  transient  techniques,  with  the  transient  techniques  yielding 
systematically  higher  values.  A  more  rigorous  analysis  of  the  transient  data  may  lead  to  additional 
details  regarding  charge  trapping/detrapping  and  relevant  timescales. 

4.3.2.  Mobility-Lifetime  Product  Measurements 

Steady-state  photoconductivity  (SS-PC)  experiments  are  commonly  applied  to  determine  the 
mobility-lifetime,  fix,  product  of  a  material.  Steady-state  photocurrent  density,  /Ph,  is  dependent 
on  photon  energy/wavelength  (/l)  and  applied  field  (£)/voltagc  (V),  as:54 

Jph&.V)  =  e0>  f  rjg(A,x)r]c(x,V)dx  [10] 

Jo 


where  e  is  elementary  charge,  &  photon  flux,  >/g  charge  carrier  generation  efficiency,  and  ;/c  charge 
collection  efficiency.  At  low  electric  fields,  photoconductivity,  crph  (Arr  =  criiiuminated  -  Ociark),  can  be 
related  directly  to  fix  via: 

Oph  =  A  a  —  (Tin  —  erdark  =  eAn/r  =  eG[ir  [11] 


Here,  the  change  in  free  charge  carrier  concentration,  An  that  occurs  upon  steady-state  illumination 
is  dependent  on  the  balance  between  the  charge  carrier  generation  rate,  G,  and  the  relaxation 
(trapping/recombination)  rate,  defined  by  the  lifetime,  r  (note  that  in  the  above  equation,  n,  /q  and 
x  represent  values  for  either  holes  or  electrons).  Generation  rate  can  be  determined  relatively 
straightforwardly  from  knowledge  of  the  photon  flux,  and  optical  response  of  the  device  via: 


G 


=  VKi-/0 


(1  —  expad) 
d 


[12] 


Where  d  represents  film  thickness,  and  a  and  R  the  absorption  coefficient  and  reflectivity  of  the 
film,  respectively.  In  this  analysis,  rjs  can  either  be  determined  independently,  assumed  to  be  1,  or 
incorporated  into  the  mobility-lifetime  term  as  >igjux. 

Another  treatment  involves  fitting  the  voltage  (V)/field  (A)-dcpendcnt  photocurrent 
(7ph)/photocurrent  density  (/Ph),  where  ultimately  the  charge  carrier  generation  rate/efficiency  term 
becomes  absorbed  into  a  saturation  photocurrent  term,  ipho,  which  can  be  treated  as  a  fitting 
parameter.  If  sufficiently  high  fields  are  probed,  and  under  the  assumption  of  uniform  light 
absorption  throughout  the  thickness  of  the  film,  JPh  can  be  fit53  as  a  function  of  V: 
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prV 

Iph (^)  —  JphO  ^2 


utV  (  Z^L 

1  — I  1  —  exp^TV 


[13] 


Note  that  if  absorption  is  assumed  to  be  primarily  at  the  surface  of  the  device,  the  Many 
equation55,56  is  used: 


700  = 


pzbV 
Jo  d2 


1  —  exp^Tbv 


1  +  —  - 
V  p 


[14] 


which  takes  into  account  a  surface  recombination  velocity  term,  sip  (here  rb  represents  the  lifetime 
within  the  bulk).  Additional  corrections  can  be  applied  if  non-uniform  electric  fields  are 
expected.57 
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Figure  28.  Schematic  of  steady-state  photoconductivity  and  dual-beam  photoconductivity 

experimental  setups. 
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Figure  29.  Photoconductivity  experimental  setup. 


An  experimental  SS-PC  apparatus  was  installed  (Figure  28  and  Figure  29)  based  on  an  OBB 
PowerArc  75  W  Xe  arc  lamp  light  source  coupled  to  a  Jobin  Yvon  (now  Horiba)  HR  640 
monochromator  to  produce  a  monochromatic  beam.  A  focused  beam  is  guided  to  a  sample  holder 
containing  the  device  under  test,  which  is  connected  to  an  HP  6632A  voltage  source  and  Keithley 
6514  electrometer.  The  experimental  setup  was  also  adapted  for  dual-beam  photoconductivity 
(DB-PC)  experiments.  Here,  the  PowerArc  Xe  lamp  is  used  as  a  monochromatic  chopped  light 
source,  and  a  second  tungsten-halogen-quartz  (THQ)  lamp  is  used  as  a  DC  bias  light  source. 

Custom  Faraday  cages  were  machined  and  assembled  for  the  experiment  [Figure  30(a)].  These 
have  an  airtight  seal  and  are  fitted  with  valves  to  allow  the  inner  space  to  be  filled  (and  pressurized) 
with  dry  nitrogen  gas  such  that  the  devices  may  be  kept  under  an  inert  atmosphere.  They  further 
shield  the  samples  from  outside  light  and  electromagnetic  interference.  A  9  mm  diameter  hole  in 
the  bottom  of  the  cage  allows  light  to  be  directed  onto  the  device.  The  opening  is  sealed  with  a  1 
mm  thick  sapphire  window,  which  allows  UV-NIR  wavelengths  to  pass  with  minimal 
reflectance/absorbance,  while  maintaining  an  airtight  space.  BNC  adapters  allow  for  connection 
to  external  electrical  circuitry. 
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Figure  30.  (a)  Faraday  cage  containing  device  under  test;  and  (b)  heterostructure  geometry  for 

photoconductivity  experiments. 
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Figure  31.  Dark  and  illuminated  current  as  a  function  of  voltage  for  a  glass/ITO/a-BAC:H  v/Cu 
heterostructure  under  ~0.2  mW/cm2  illumination  at  450  nm. 


Heterostructures  for  measurement  were  fabricated  by  depositing  a-BAC:Hy  onto  a  glass/indium 
tin  oxide  (ITO)  substrate,  and  then  depositing  a  top  Cu  contact  by  sputtering  [Figure  30(b)],  Films 
used  in  the  results  shown  here  were  grown  at  a  temperature  of  300  °C  10  W,  0.4  Torr,  200  seem 
total  flow,  and  a  partial  Ar  +  orthocarborane  pressure  of  1  for  30  min,  resulting  in  films  with  a 
thickness  of  750  nm.  The  devices  were  then  mounted  to  Cu-coated  FR  board  with  silver  paint/gold 
wire,  and  soldered  to  the  BNC  connections  within  the  Faraday  cage.  Using  the  above-described 
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setup,  dark  and  illuminated  current  were  measured  as  a  function  of  wavelength  from  300  to  800 
nm  and  of  voltage  from  0  to  5  V.  Photocurrent  was  taken  as  the  difference  between  the  two,  from 
which  was  extracted  photoconductivity  via: 


dl 

°vh  =  \ Ta 


[15] 


Generation  rate  was  calculated  via  Equation  12.  Photon  flux  was  calculated  from  a 
measurement  of  power  output  at  the  sample  position  (after  the  Faraday  cage  to  take  into  account 
losses  due  to  the  sapphire  window),  which  was  normalized  to  irradiation  spot  size.  Quantum 
efficiency  was  assumed  to  be  1.  Absorption  coefficient,  a,  was  determined  from  UV-Vis 
absorbance  data  [Figure  32(a)],  and  reflectivity,  R,  assumed  to  be  0.  Note  that  both  the  quantum 
efficiency  and  reflectivity  assumptions  are  oversimplifications  of  the  analysis.  Mobility-lifetime 
product  calculated  based  on  this  generation  rate  using  Equation  11  varies  as  a  function  of 
wavelength,  peaking  at  ~2.8  eV  (450  nm),  at  which  point  it  was  calculated  to  be  2  x  10-14  cm2/V. 
A  qualitative  snapshot  of  photoconductivity  as  a  function  of  illumination  energy  is  provided  by 
the  aph/0  curve  [Figure  32(b)],  wherein  it  can  be  seen  that  it  roughly  mirrors  the  absorption  curve, 
though  the  photoconductivity  onset  is  at  about  1  eV  lower  energy.  SS-PC  measurements  on  a 
number  of  samples  grown  under  similar  conditions  yielded  consistent  jut  values  on  the  order  of  10" 
13 — 10_14  cnr/V.  Mobility  measurements  on  similar  films  yielded  n  of  the  same  order  of  magnitude, 
suggesting  lifetimes  on  the  order  of  1  s. 


Figure  32.  (a)  Absorption  coefficient  for  a-BYC:Hv  on  glass/ITO.  (b)  Photoconductivity 
(normalized  to  photon  flux)  as  a  function  of  energy. 
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4.4.  Electronic  Structure  and  Charge  Transport  Models 

4.4.1.  Electronic  Structure  Model 
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Figure  33.  Schematic  for  the  density  of  states,  band  diagram,  and  absorption  spectrum  for  an 

amorphous  solid. 


Toward  understanding  the  charge  transport  properties  of  a  material,  it  is  instructive  to  analyze 
its  electronic  structure  (Figure  33).  Amorphous  semiconductors  are  unique  from  crystalline 
semiconductors  in  that  they  exhibit  band  tails  arising  from  disorder-induced  localized  states.  These 
exponential  band  tails,  known  as  Urbach  tails,  are  separated  from  the  conduction  and  valence  bands 
by  the  mobility  edges,  which  define  the  mobility  gap — analogous  to  the  band  gap  in  crystalline 
solids.  While  the  localized  states  that  comprise  the  Urbach  tails  can  be  thought  of  as  a  set  of 
shallow  traps,  amorphous  solids  may  also  contain  deep  traps  or  mid-gap  states.  A  powerful 
technique  for  investigating  electronic  structure  is  optical  absorption  spectroscopy,  where  the 
absorbance  spectrum  represents  a  superposition  of  optical  transitions  between  combinations  of 
extended,  localized,  and  mid-gap  states.58,59  In  a  plot  of  absorption  coefficient  vs  energy,  three 
main  regions  can  be  identified:  region  A,  at  higher  energies  and  with  higher  absorption  intensities, 
is  associated  with  transitions  between  extended  states;  region  B,  at  moderate  energies  and 
absorption  intensities,  is  associated  with  transitions  to  or  from  localized  tail  states,  and  has  been 
shown  to  be  primarily  defined  by  transitions  between  localized  valence  band  states  and  extended 
conduction  band  states;60  and  region  C,  at  the  lowest  energy  and  absorption  intensities,  is 
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associated  with  transitions  involving  mid-gap  defect  states.  Tauc  showed  that  region  A  can  be 
modeled  by  a  power  law  according  to: 


B(E  —  Eg)n 


[16] 


which  can  be  rearranged  to  (with  n  defined  as  2,  as  is  typical  for  most  amorphous  solids): 

(ocE)i  =  b\{E  -  ETauc)  [17] 

such  that  the  Tauc  optical  band  gap,  Txauc,  and  the  Tauc  parameter,  Bu2,  can  be  extracted  from  a 
plot  of  ( aE)u 2  vs  E  as  the  intercept  and  slope  of  the  linear  region,  respectively.  Other  methods  have 
been  employed  to  determine  the  optical  band  gap  of  amorphous  solids,61  including  the 
isoabsorption  method,60-62  where  a  gap  is  empirically  defined  at  a  given  absorption  coefficient 
value  (typically  Em  at  1  x  104  cm-1),  as  well  as  the  use  of  alternative  values  for  “n”  in  the  Tauc 
expression  (typically  3,  as  may  be  appropriate  for  boron-rich  solids).63,64  Region  B  can  be  fit  to  the 
expression: 

a  —  a0  exp  )  [18] 


where  the  Urbach  energy,  Eu,  can  be  obtained  as  the  slope  of  the  curve  for  a  plot  of  ln(a)  vs  ln(£). 


An  analysis  of  the  absorbance  spectra  for  a  large  number  of  a-BAC:Hv  films  was  completed. 
Details  of  the  analysis  are  described  in  Nordell  et  al.  JAP  118  (2015)  035703,  and  a  subset  of  the 
results  can  be  found  in  this  publication  as  well  as  a  subsequent  one  Nordell  et  al.  (2015) 
submitted,  MATCHEMPHYS-D-15-01944.  Three  different  methods  were  applied  to  extract  optical 
band  gap:  a  Tauc  treatment  with  n  =  2  and  n  =  3  as  well  as  an  isoabsorption  treatment  with  Em  defined 
as  the  energy  at  which  a  -  1  x  104  cm-1.  It  was  found  that  isxauc  (n  =  2)  and  Em  values  were  relatively 
close,  and  could  both  serve  as  appropriate  measures  of  optical  band  gap  [see  Figure  34(a)  for  a 
comparison  of  the  three].  The  three  electronic  structure  parameters,  E% ,  BU2,  and  Eu  were  found  to 
range  widely  as  a  function  of  thin-film  preparation  and  composition,  with  optical  band  gap  ranging 
from  1. 7-4.1  eV,  B112  from  100-2200  eV-1/2  cm-1/2,  and  E\\  from  100-800  meV. 
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Figure  34.  Correlation  between  (a)  band  gap  values,  £Vauc  [n  =  2  (filled  circle)  and  n  =  3  (blue 
triangle]  and  Eo4  (open  circle),  and  Urbach  energy  and  (b)  Tauc  slope  parameter,  BU2,  and 
Urbach  energy,  for  a  series  of  a-BC:H  films  [Figure  from  Nordell  et  al.  (2015)  submitted, 

MATCHEMPHYS-D-15-019441 . 


Band  gap  is  correlated  to  hydrogen  concentration  (Figure  2  and  Figure  3),  where  a  hydrogen- 
rich  low-density  film  has  a  high  band  gap  and  vice  versa.  This  is  the  same  behavior  observed  in  a- 
Si:H,  one  explanation  for  which  pertains  to  the  decrease  in  network  connectivity  associated  with 
the  incorporation  of  one-fold-coordinate  hydrogen  atoms.  The  Urbach  energy  and  Tauc  parameter, 
both  measures  of  disorder  in  an  amorphous  solid,  are  also  correlated  with  hydrogen  concentration: 
Urbach  energy  decreases  with  increasing  hydrogen  concentration  ( Bu 2  shows  the  opposite  trend, 
but  for  this  parameter  high  values  indicate  low  disorder),  which  is  the  same  behavior  observed  in 
a-Si:H,  and  can  be  explained  by  a  decrease  in  network  strain  associated  with  hydrogen 
incorporation,  and  thus  decrease  in  structural  disorder.59'65  We  also  found  in  Nordell  et  al.  (2015) 
submitted,  MATCHEMPHYS-D-15-01944  that  the  carbon  content  plays  a  smaller  role  in  defining 
electronic  structure,  even  though  the  major  trends  can  be  associated  with  hydrogen  content. 
Finally,  we  see  that  band  gap  and  Urbach  energy  are  inversely  correlated  to  one  another  in  a  linear 
fashion  [Figure  34(a)],  with  larger  band  gaps  associated  with  smaller  band  tails,  and  vice  versa. 
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Overall,  we  have  found  a  remarkable  correlation  between  electronic  structure  parameters  and 
chemical  composition  over  a  wide  range  of  stoichiometries.  Although  relatively  simple,  this 
electronic  structure  model  provides  a  framework  for  interpreting  gross  trends  in  charge  transport 
in  order  to  zero  in  on  the  regime  appropriate  for  fabricating  a-BYC:Hv  detectors. 


4.4.2.  Charge  Transport  Model 
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Figure  35.  Correlations  between  (a)  mobility  and  Tauc  band  gap  (n  =  2),  (b)  mobility  and 
Urbach  energy,  (c)  Tauc  band  gap  and  Urbach  energy,  and  (d)  mobility  and  Es/E\j. 


The  charge  transport  parameters  were  compared  against  electronic  structure  parameters  to 
investigate  correlations.  These  correlations  were  discussed  in  some  detail  for  a  subset  of  films  in 
Nordell  et  al.  (2015)  submitted,  MATCHEMPHYS-D-15-01944.  What  we  observed  is  that  the 
charge  transport  parameters  do  not  show  smooth  correlations  with  electronic  structure  and 
chemical  composition  properties,  but  rather  appear  to  exhibit  different  correlated  regimes.  In 
Figure  35(a),  we  observe  that  mobility  exhibits  a  relatively  strong  inverse  linear  correlation  with 
band  gap  below  ~3  eV,  and  what  appears  to  be  a  weaker  direct  linear  correlation  above  this  value. 
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The  same  type  of  behavior  is  observed  in  the  mobility  vs  Urbach  energy  curve,  with  a  threshold 
of  -500  meV  [Figure  35(b)].  We  have  already  discussed  how  Urbach  energy  and  band  gap  are 
inversely  linearly  correlated.  This  is  illustrated  in  Figure  35(c)  for  a  much  larger  set  of  samples. 
Figure  35(d)  demonstrates  how  juf  can  be  correlated  to  both  parameters  via  a  joint  term  EJEy  : — in 
the  low  iig/high  Ew  region,  the  correlation  is  particularly  clear. 

Unlike  other  thin-film  properties  explored,  which  tend  to  exhibit  linear  correlations  with  each 
other  and  with  chemical  composition,  the  unique  trends  demonstrated  by  the  charge  transport 
properties  may  suggest  the  existence  of  more  than  one  relevant  transport  mechanism.  We  initially 
hypothesized  that  mobility  would  trend  with  Urbach  energy,  specifically  decreasing  with  increasing 
E\j  and  thereby  increasing  film  disorder,  as  in  the  case  of  amorphous  hydrogenated  silicon.  However, 
we  find  that  maximum  mobility  values  are  found  at  maximum  E\j  values,  and  that,  at  low  band  gap, 
mobility  decreases  with  decreasing  E\j.  Two  different  charge  transport  mechanisms  can  be  considered 
for  an  amorphous  semiconductor.  In  a  multiple  trapping  (MT)  and  release  mechanism,  charge  transport 
occurs  via  extended  band  states,  and  mobility  is  reduced  through  trapping  by  localized  states.  In  this 
scenario,  a  higher  Urbach  energy  would  imply  wider  band  tails,  greater  trapping,  and  thus  lower 
mobility.  In  variable  range  hopping  (VRH)  conduction,  charge  transport  occurs  due  to  transitions 
between  localized  states  within  the  gap,  and  conductivity  is  thus  related  to  the  transition  probability 
for  hopping  between  sites,  and  can  ultimately  be  related  to  the  concentration  of  sites. 

Figure  36  distinguishes  three  different  electronic  structure  regimes.  Films  with  high  band  gaps 
also  have  low  Urbach  energies.  This  regime  is  associated  with  lower-density,  hydrogen-rich  films, 
which  should  have  low  strain  and  low  disorder — hence  the  narrow  band  tails.  In  this  regime, 
decreases  in  Eg  and  increases  in  E\j  are  associated  with  a  trend — albeit  scattered — toward  lower 
mobilities.  This  trend  would  be  consistent  with  a  multiple  trapping  mechanism,  where  a  greater 
number  of  localized  states  leads  to  increased  trapping  and  therefore  decreased  mobility.  Films  with 
low  band  gaps  have  very  high  Urbach  energies.  In  fact,  the  band  tail  width  relative  to  band  gap 
implies  that  the  exponential  distribution  of  localized  states  fills  the  entire  gap.  This  regime  is 
associated  with  high-density,  hydrogen-deficient  films.  From  rigidity  theory  arguments,65  these 
amorphous  films  would  be  expected  to  be  significantly  overconstrained,  and  thus  highly  strained, 
which  would  explain  the  very  high  disorder  parameters.  In  this  electronic  structure  regime, 
increases  in  Eg  and  decreases  in  E\j  are  associated  with  a  trend  toward  lower  mobilities.  One 
explanation  for  this  would  be  that  a  variable  range  hopping  transport  mechanism  dominates,  where 
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the  main  conduction  pathway  is  via  localized  rather  than  band  states.  This  is  reasonable  due  to  the 
high  density  of  localized  states  within  the  gap.  In  this  case,  mobility  would  increase  with  a  greater 
density  of  localized  states,  as  opposed  to  in  a  MT  scenario.  These  transport  models  are  still  very 
crude,  but  they  provide  an  initial  framework  for  thinking  about  the  structure-property  relationships 
in  a-B.vC:Hv  films.  For  example,  a  goal  to  “increase  mobility  by  decreasing  Urbach  energy,”  as 
would  be  the  case  in  a-Si:H,  would  not  hold  in  higher  density  a-BAC:Hv  where  mobility — and 
Urbach  energy — are  highest.  In  this  region,  it  may  be  important  to  attempt  to  distinguish  types  of 
traps  (e.g.,  shallow  vs  deep)  or  types  of  disorder  (e.g.,  structural  vs  substitutional)  in  thinking  about 
further  optimizing  charge  transport  properties.  Overall,  more  studies  will  be  needed  to  confirm 
hypothesized  models  (e.g.  variable  temperature  conductivity  measurements)  and  trapping  profiles. 


F  (eV) 


Multiple  trapping 
Higher  mobility 


MT/VRH 
Low  mobility 


Variable  range  hopping 
Higher  mobility 

Vanable  range  hopping  (high  T) 


(c)  Spatial  coordinate 


Figure  36.  Three  different  proposed  electronic  structure  models  and  associated  charge  transport 
mechanisms  in  a-BYC:Hy:  (A)  high  band  gap  and  low  Urbach  energy  associated  with  a  multiple 
trapping  mechanism,  (C)  low  band  gap  and  high  Urbach  energy  associated  with  a  variable  range 
hopping  mechanism,  and  (B)  an  intermediate  regime. 
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4.5.  Detector  Modeling 

4.5.1.  Theory 

Our  goal  was  to  model  the  detection  efficiency,  specifically  for  the  case  of  a  thin-film  detector 
with  relatively  low  mobility,  in  order  to  determine:  (a)  what  detector  specifications  and  a-BxC:H y 
figures  of  merit  are  needed  to  achieve  a  targeted  efficiency,  and/or  (b)  given  known  detector 
specifications  and  a-BvC:Hv  figures  of  merit,  what  maximum  detection  efficiency  is  theoretically 
possible.  Because  detectors  are  traditionally  based  on  greater  thicknesses  of  bulk  material  with 
higher  mobilities,  the  associated  theory  and  assumptions  cannot  necessarily  be  applied  to  the  type 
of  device  proposed  here  and  must  be  reviewed. 


Thickness  (jam) 

Figure  37.  Absorption  probability  of  thermal  neutrons  as  a  function  of  B-tC  thickness. 


The  first  step  in  determining  detection  efficiency  involves  understanding  the  interaction  of 
neutrons  with  a  material.  The  macroscopic  cross  section,  defined  as  the  interaction  probability  per 
unit  length  of  an  incident  particle  in  a  medium,  can  be  calculated  for  any  composition  of  matter 
based  on  the  microscopic  cross  sections  of  the  atomic  constituents.  For  a  boron-based  detector,  we 
are  interested  in  the  interaction  between  10B  and  thermal  neutrons  (i.e.,  those  with  an  energy  of 
-0.025  eV,  the  energy  corresponding  to  room  temperature),  where  10B  has  a  very  high  microscopic 
thermal  neutron  cross  section  of  3840  bam66  (1  barn  =  10-24  cm2).  If  we  assume  a  material  with  a 
stoichiometry  of  B4C  and  density  of  2.52  g/cm3,67  the  macroscopic  cross-section  for  a  naturally 
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abundant  sample  (-20%  of  10B)  is  84.38  cm-1  and  for  an  enriched  sample  (100%  10B)  is  421.92 
cm'1.  The  corresponding  penetration  depths  (i.e.,  the  reciprocal  of  the  macroscopic  cross  section) 
are  118.5  pm  and  23.7  pm,  respectively.  The  absorption  probability  for  thermal  neutrons  as  a 
function  of  B4C  thickness  are  given  in  Figure  37.  It  is  evident  that  for  an  enriched  sample,  nearly 
maximal  neutron  absorption  can  be  achieved  for  thin  100  pm  films,  and  only  30  pm  of  material  is 
sufficient  for  -75%  absorption. 

After  neutron  absorption,  the  next  phenomenon  that  must  be  considered  is  the  fate  of  the 
neutron  capture  reaction  charged  particle  products.  The  absorption  of  a  neutron  by  a  10B  nucleus 
results  in  the  emission  of  a  positively  charged  alpha  particle  and  a  negatively  charged  lithium  ion. 
About  94%  of  this  interaction  leaves  the  Li  ion  in  the  excited  state  with  the  energy  distribution  Eu 
=  0.84  MeV  and  Ea  =  1.47  MeV,  and  the  remaining  6%  interaction  leaves  the  Li  ion  in  the  ground 
state  with  the  energy  distribution  Eu  =  1.02  MeV  and  Ea  =  1.78  MeV.66  Because  of  the  negligible 
momentum  of  the  incident  thermal  neutron  relative  to  that  of  the  emitted  particles,  the  emitted 
particles  move  opposite  to  one  another  with  a  random  emission  angle.  As  long  as  the  emitted 
particles  move  within  the  detector  material  (B4C),  they  excite  electron-hole  pairs  along  their  paths. 
If  U  is  the  total  energy  deposited  by  the  emitted  particles  and  Ev  is  the  energy  required  to  create  an 
electron-hole  pair,  the  number  of  electron-hole  pairs  (Ne hP)  generated  will  be: 

Nehp  =  [19] 

bP 

If  an  electric  field  E  is  applied  across  a  detector  material  in  a  capacitor  geometry,  the  excited  e-h 
pairs  are  displaced  against  the  field  and  induce  a  current  in  the  external  measurement  circuit  given 

as:68’69 


/  = 


qq.E 


[20] 


where  q  =  Ue/Epis  the  charge  deposited  (charge  of  excited  electrons/holes),  fx  is  the  charge  carrier 
mobility  in  the  detector  material,  and  d  is  the  detector  thickness  (i.e.,  electrode  spacing).  The 
current  can  then  be  integrated  for  a  certain  time,  t,  with  the  total  charge  collection  represented  by: 


Qcoll 


1  utE 

Idt  —  q - 


d 


[21] 
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In  a  real  detector,  the  excited  charge  carriers  are  distributed  along  the  detector  thickness 
according  to  the  energy  deposited  by  the  emitted  Li  ion  and  a  particle.  The  number  of  excited 
carriers  decreases  exponentially  with  time  as  a  result  of  carrier  trapping  and/or  recombination,  and 
the  carrier  concentration  at  any  time  t  is  given  as: 

n(x  +  vt,  t)  =  n(x,  0)  exp  ^ j  [22] 

where  n(x,  0 )  is  the  initial  excited  carrier  concentration  at  x  and  x  is  the  carrier  lifetime.  Because 
of  the  decreasing  carrier  concentration,  the  induced  current  also  decays  exponentially  and  the 
collected  charge  corresponding  to  the  charge  deposited  in  a  small  thickness  dx  at  x  within  the 
integration  time  t  becomes: 

dQcon  =  f  q(x,0)  dx  ^-exp(-^  dt  [23] 

where  q(x,  0)dx  is  the  charge  deposited  within  a  small  thickness  dx  and  is  obtained  from  the 
stopping  power,  S(x),  of  the  detector  material  for  a  particles  and  Li  ions  as: 

Six') 

q(x,0)  =  -—e  [24] 

hp 


The  total  charge  collected  is  then  obtained  as: 


Qcoll 


'  d  L 


q(x,  0)  dx 


IexpK) 


dt 


[25] 


where  xmax  -  Xmin  is  the  range  within  which  the  charge  carriers  are  excited  by  the  neutron  absorption 
event. 


To  collect  the  maximum  deposited  charge,  the  minimum  integration  time  must  be  equal  to  or 
greater  than  the  carrier  transit  time,  ttI  =  dl/iE.  Detectors  based  on  higher  mobility  materials  (e.g., 
Si)  can  have  low  integration  times  (e.g.,  on  the  order  of  microseconds)  because  of  the  short  transit 
time  of  the  carriers  even  in  the  case  of  thick  detectors  (e.g.,  on  the  order  of  ones  of  millimeters). 
For  low  mobility  materials,  the  transit  time  becomes  long,  which  sets  a  higher  minimum  on  the 
required  integration  time.  However,  the  integration  time  cannot  always  be  increased  beyond  a 
certain  limit.  This  limit  might  be  set  by  the  time  gap  between  successive  detection  events — that  is, 
if  the  charge  is  integrated  over  too  long  of  a  period,  two  or  more  detection  events  may  result  in  an 
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overlapped  signal  pulse  (note  that  this  may  not  be  a  problem  in  all  detection  scenarios,  but  to  obtain 
a  neutron  detection  efficiency,  individual  neutrons  must  be  counted).  For  spectroscopic  studies, 
the  match  up  of  integration  and  transit  time  is  also  important.  Additionally,  from  a  practical 
experimental  standpoint,  the  integration  time  may  be  limited  by  the  electronics  available,  which 
are  typically  tailored  to  higher  mobility  materials  and  shorter  timescale.  Thus,  because  it  may  be 
necessary  or  desirable  to  use  an  integration  time  smaller  than  the  transit  time,  one  must  take  into 
account  the  incomplete  charge  collection  that  results  from  this  scenario.  If  the  integration  time 
used  is  less  than  the  transit  time,  a  boundary  is  created  at: 

x0  —  d  —  qtE  [26] 

where  all  the  charge  carriers  generated  at  x  >  xo  can  reach  the  collecting  electrode  within  the 
integration  time  t,  whereas  those  generated  at  x  <  xo  can  cover  a  maximum  distance  of  filE.  Thus 
the  distance  that  a  charge  carrier  travels  becomes  smaller  than  the  detector  thickness  in  both  of 
these  cases  and  leads  to  incomplete  charge  collection.  The  total  charge  collected  in  a  typical  case 
of  Xmin  <  xo  <  Xmax  can  therefore  be  written  as  the  sum  of  two  terms  as: 


LITE  r  /  t\  Cx°  U.tE  fXmax  r  /  d  —  x\i 

<2coii  =  -y  [1  -  exp  (- -JJ  I  q(x)dx  +  —\  q(x)dx  ^1  -  exp  ^ - — jj  [27] 

Jxm\n  JX0  A- 


Figure  38.  Schematic  of  neutron  detection  in  a  thin-film  device. 
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On  the  basis  of  the  above  theory,  a  numerical  calculation  scheme  for  the  detection  efficiency 
of  an  a-BxC:Hy  detector  was  implemented  using  Python.  In  setting  up  the  detection  scenario,  we 
assume  that  neutrons  are  incident  normally  on  the  surface  of  a  planar  thin-film  detector  of  thickness 
cl,  and  a  neutron  is  absorbed  by  a  10B  nucleus  at  a  distance  x  from  the  incident  surface  with  a 
probability  density  P(x )  determined  by  the  macroscopic  cross-section  (Figure  38).  If  an  alpha 
particle  is  emitted  at  an  angle  6  with  respect  to  the  field  direction,  a  Li  ion  is  emitted  at  an  angle  jt 
+  6.  The  charge  carriers  are  thus  generated  between  xmin  andxmax  determined  by  the  emission  angle 
6  and  the  ranges  Ra  and  Ru  of  the  a  particle  and  Li  ion,  respectively.  If  this  range  lies  inside  the 
detector,  the  total  energy  of  the  reaction  products  is  deposited  and  leads  to  the  excitation  of  charge 
carriers;  however,  if  the  reaction  products  escape  the  detector,  only  a  fraction  of  the  total  energy 
is  deposited. 

4.5.2.  Calculation  and  Assumptions 

The  detection  efficiency  calculation  comprises  the  following  general  steps: 

1.  A  random  number  x  is  generated  with  probability  density: 

P(x)  =  E  exp(— Ex)  ;  x  >  0  [28] 

where  E  is  the  macroscopic  cross-section.  The  random  number  x  represents  the  interaction 
position  of  an  incident  neutron  with  10B  inside  the  detector.  Any  values  with  x  >  cl  are  discarded 
and  the  event  is  counted  as  not  detected. 

2.  Another  random  number  6  c  [0,  71]  is  generated,  which  represents  the  emission  angle  of  the 
alpha  particle  from  the  interaction  position  x. 

3 .  The  range  xmin  -  Xmax  over  which  the  energy  of  the  a  particle  and  Li  ion  is  deposited  is  calculated 
based  on  the  interaction  position  x  and  emission  angle  6. 

4.  The  boundary  xo  is  calculated  based  on  the  defined  values  of  mobility,  /i,  integration  time,  t, 
and  applied  electric  field,  E. 

5.  The  charge  deposited  per  unit  distance,  q(x),  is  calculated  based  on  the  known  stopping  power 
of  the  detector  material  and  the  emission  angle  9. 

6.  The  charge  collected  by  the  external  circuit  is  calculated  using  the  following  expressions: 


For  Xmin  >  Xo: 
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Qcoll  — 


/ttE 


f 

J  X. 


g(x)dx 


[29] 


For  Xmax  <  X( !  ■ 


Qcoll  -  Q 


[30] 


For  Xmin  <  X0  <  Xmax- 


Qcoll 


[irE 


/  t\]  fX°  [IT E  fXmax 

1  —  exp  I  —  1  I  q(x)dx  +  -j-  I  q(x)dx 

-*-min  *0 


[31] 


Here,  for  a  hole-only  device  where  the  mobile  holes  are  moving  along  the  x-axis  (i.e.,  along 
the  electric  field  axis),  hr  =  (d  -  x)//uE  and  Q  is  the  total  charge  deposited. 

7.  The  collected  charge  is  compared  against  the  minimum  detection  threshold:  if  it  exceeds  the 
threshold,  the  event  is  counted  as  detected. 

8.  The  total  detection  efficiency  is  then  calculated  as: 

Number  of  detected  counts 

Detection  Efficiency  = - - - - - - - - - x  100%  [32] 

Total  number  of  incident  neutrons 


For  the  calculation  of  the  detection  efficiency  in  an  a-BvC:Hv  detector,  the  following  specific 
assumptions  were  applied.  We  note  that  these  assumptions  are  expected  to  provide  a  close 
approximation  to  the  true  behavior  of  an  a-BAC:Hv  device,  but  can  be  modified  in  a  plug-and-play 
fashion  to  reflect  new  information  regarding  the  properties  of  the  material/device  under  test  and/or 
the  specific  experimental  setup  and  limitations. 

1.  The  detector  material  was  assumed  to  have  a  composition  of  B4C  and  density  of  2.52  g/cm3. 
In  reality,  the  material  used  will  have  a  B/C  ratio  between  3.5  and  5,  5-20%  hydrogen 
concentration,  and  a  density  between  1.8  and  2.4  g/cm3.  Future  treatments  can  be  modified  to 
include  the  precise  composition  of  the  bulk  detector  material. 
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2.  The  energy  deposition  for  an  a  particle  and  a  Li  ion  in  the  B4C  material  was  obtained  using  the 
SRIM  (Stopping  and  Range  of  Ions  in  Material)70  simulation  package.  The  stopping  powers 
(i.e.,  the  energy  deposited  per  unit  distance)  for  the  particles  are  shown  in  Figure  39. 


Figure  39.  Bragg  curves  for  a  1.47  MeV  a  particle  and  a  0.84  MeV  Li  ion  in  B4C. 


From  the  SRIM  calculation,  the  range  for  an  a-particle  in  B4C  was  found  to  be  3.32  pm  and 
for  a  Li  ion  to  be  1.69  pm.  For  simplicity,  the  stopping  power  of  the  a  particle  was  approximated 
to  be  uniform  over  its  range  by: 

keV 

Sa  (x)  =  440.12  -  [331 

pm  L  J 

and  that  of  the  Li  ion  was  approximated  by  a  linear  equation  as: 

keV 

SLiO)  =  -433.2x  +  884  —  [34] 

3.  The  average  energy  required  to  produce  one  e-h  pair  is  given  by:71 

14 

Ep  —  —Eg  +  r  )  0.5  <  r  <  1.0  eV  [35] 

where  Eg  is  the  band  gap  of  the  material.  Based  on  a  representative  band  gap  value  of  3-3.5 
eV,  the  energy  cost  per  e-h  pair  was  approximated  to  be  10  eV. 

4.  Only  one  type  of  charge  carrier — specifically,  holes — was  assumed  to  be  mobile,  and  the 
direction  of  the  applied  electric  field,  and  thus  carrier  drift,  are  in  the  direction  of  the  incident 
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neutrons.  The  choice  of  holes  as  mobile  carriers  is  supported  by  evidence  of  p-type 
conductivity  in  crystalline72  and  amorphous22  boron  carbide,  and  has  been  consistent  with 
charge  transport  measurements  within  our  group.  This  should  represent  a  lower  minimum  on 
collection  efficiency  for  the  thin-film  case. 

5.  The  minimum  charge  required  for  detection  is  set  at  10  fC  within  an  integration  time  of  1  ms. 
This  leads  to  an  average  minimum  current  of  10  pA  through  the  detector  which  should  be 
detectable  with  sensitive  electronics. 

6.  The  range  of  parameters  used  for  the  calculations  were: 

•  Drift  velocity  (juE):  5  x  10“3  -  50  cm/sec  with  20  data  points  per  decade 

•  Lifetime  (r):  10'5-1  sec  with  1  data  point  per  decade 

•  Thickness  ( d ):  1-20  pm  in  equal  intervals  of  1  pm 

7.  The  collected  charge  for  the  long  lifetime  case  depends  on  the  integration  time,  that  is: 

[itE 

Qcoll  —  ^  ^deposited  [36] 

In  order  to  truncate  the  calculation,  combinations  of  fiE,  t,  and  d,  for  which  ulE/d  <0.1  were 
not  used  as  these  combinations  yield  a  charge  collection  value  well  below  the  threshold  for 
detection.  Similarly  for  the  short  lifetime  case,  the  combinations  for  which  firEld  <0.1  were 
also  discarded. 

8.  The  calculation  is  based  on  10,000  neutrons  incident  on  the  detector,  and  the  area  of  the 
detector  is  assumed  to  be  very  large. 

9.  The  charge  collection  efficiency  was  calculated  based  on  the  total  energy  collected  relative  to 
the  total  energy  deposited  (for  all  10,000  neutrons),  however  the  detection  efficiency  was 
calculated  independently  for  each  individual  case. 

4.5.3.  Results  and  Discussion 

Figure  40  shows  the  charge  collection  efficiency  as  a  function  of  carrier  drift  velocity  (juE)  for 
different  detector  thicknesses  and  carrier  lifetimes.  For  short  lifetimes,  the  charge  collection 
efficiency  increases  slowly  and  saturates  with  increasing  drift  velocity  because  charge  carrier 
trapping  on  timescales  faster  than  the  integration  time  leads  to  incomplete  charge  collection.  For 
longer  lifetimes,  the  charge  collection  efficiency  saturates  at  a  higher  rate  because  most  of  the 
deposited  charges  are  collected  before  they  are  trapped.  Further,  it  is  evident  that  the  drift  velocity 
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must  be  increased  in  proportion  with  the  thickness  to  obtain  equivalent  charge  collection 
efficiency,  that  is,  either  mobility  or  electric  field  needs  to  increase. 
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Figure  40.  Charge  collection  efficiency  (CCE)  as  a  function  of  drift  velocity  (v  =  juE)  for  a  range 
of  a-BvC:Hy  thicknesses  at  different  lifetimes:  (a)  10  ps,  (b)  1  ms,  and  (c)  1  s. 


It  is  also  evident  from  Figure  40  that  the  charge  collection  efficiency  saturates  at  50%  for  thin 
films,  and  gradually  increases  above  this  value  with  increasing  film  thickness.  For  thin  films,  the 
charge  collection  efficiency  determined  here  is  very  closely  approximated  by  the  charge  collection 
efficiency  for  the  uniform  absorption  treatment  in  steady-state  photocurrent  theory: 


utE  I  utE 

CCE  =  !-r- 1  1  -  '—r-  (  1  -  exp 


inEj 


d 


d 


[37] 
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where  the  maximum  theoretical  efficiency  is  limited  to  50%  for  high  n rE/d  (because  we  have 
assumed  single-carrier  transport).  In  the  uniform  absorption  case,  those  carriers  generated  near  the 
incident  surface  and  drifting  the  full  length  of  the  detector  will  result  in  100%  charge  collection, 
whereas  those  generated  at  the  far  surface  with  negligible  distance  to  travel  before  reaching  the 
collection  electrode  will  result  in  0%  charge  collection — thus,  giving  an  average  of  50%  collection. 
When  the  thickness  of  the  detector  medium  becomes  very  large  compared  to  the  penetration  depth 
of  the  incident  particles,  the  surface  approximation  treatment  may  be  applied,  which  predicts  a 
charge  collection  efficiency  of  100%  (i.e.,  all  carriers  are  generated  at  the  incident  surface,  and 
drift  the  full  length  of  the  detector  to  yield  100%  charge  collection).  In  the  detector  simulations 
presented  here,  the  system  moves  away  from  a  purely  uniform  absorption  scenario  towards  a 
surface  absorption  scenario,  but  even  for  the  thickest  films  modeled  (20  pm),  the  maximum  CCE 
remains  much  lower  than  100%. 


Figure  41.  Charge  collection  efficiency  (CCE)  as  a  function  of  mean  carrier  drift  length  (2  = 
inE)  for  a  range  of  a-BvC:Hv  thicknesses  at  different  lifetimes:  (a)  10  ps,  (b)  1  ms,  and  (c)  1  s. 
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Figure  41  shows  the  variation  of  charge  collection  efficiency  as  a  function  of  mean  carrier 
drift  length  (/.  =  jutE)  for  different  lifetimes.  These  graphs  are  very  similar  to  those  in  Figure  40, 
but  here  it  is  made  clear  that  a  higher  value  of  mobility-lifetime  product  (jut)  is  not  necessarily 
sufficient  for  maximum  charge  collection  in  all  detection  scenarios,  although  it  is  generally  the 
case  for  conventional  detectors  with  high  mobility  materials  requiring  short  integration  times. 
Here,  for  short  lifetimes  [Figure  41(a)  and  (b)],  there  is  little  variation  in  charge  collection 
efficiency  with  drift  length.  However  in  the  case  of  long  lifetimes  [Figure  41(c)],  due  to  a 
corresponding  low  mobility,  the  actual  distance  travelled  by  the  carriers  within  the  integration  time 
is  in  fact  less  than  the  drift  length,  which  makes  the  charge  collection  efficiency  low  even  for  large 
values  of  L  In  other  words,  drift  length — the  average  distance  that  the  carriers  can  travel  if  they 
are  allowed  sufficient  time — does  not  necessarily  represent  the  true  distance  that  carriers  travel 
within  a  specified  integration  time  less  than  their  lifetime,  so  higher  drift  length  (i.e.,  higher  jut 
product)  does  not  guarantee  high  charge  collection  and  one  must  also  consider  the  values  of  ji  and 
r  independently. 

Because  we  are  interested  in  detection  efficiency,  achieving  100%  charge  collection  efficiency 
is  not  necessarily  the  goal,  but  rather  we  seek  to  achieve  a  CCE  sufficient  to  meet  the  minimum 
threshold  for  registering  a  detection  event.  In  Figure  42,  detection  efficiency  is  plotted  as  a  function 
of  CCE  for  a  range  of  detector  thicknesses  at  different  carrier  lifetimes.  A  common  feature  in  the 
plots  is  the  onset  of  a  detection  threshold  at  a  little  over  20%  CCE.  Once  the  threshold  is  crossed, 
the  detection  efficiency  increases  slowly  to  the  saturation  value  for  thin  detectors  and  short  carrier 
lifetimes  and  rises  quickly  to  the  saturation  value  for  thick  detectors  and  long  carrier  lifetimes. 
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Figure  42.  Detection  efficiency  as  a  function  of  charge  collection  efficiency  for  a  range  of 
detector  thicknesses  at  different  lifetimes:  (a)  10  ps,  (b)  1  ms,  and  (c)  1  s. 
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Figure  43.  Detection  efficiency  as  a  function  of  fix  Eld  and  thickness  at  different  lifetimes:  (a)  x  < 

1  ms,  (b)  x  =  10  ms,  (c)  x  =100  ms,  and  (d)  x  =  1  s. 


Figure  43  shows  the  detection  efficiency  as  a  function  of  /y xE/d  and  thickness  for  different 
carrier  lifetimes.  For  small  lifetimes  (r  <  t),  the  detection  onset  occurs  when  /ixE/d  ~  0.35.  For 
lifetimes  greater  than  the  integration  time  (r  >  t ),  the  detection  onset  occurs  when  uE/d  ~  350, 
which  yields  a  current  of  ~13  pA  for  the  maximum  possible  charge  collected.  The  approximate 
mobility  required  to  achieve  neutron  detection  is  then: 

d 

[x  >  350  x  —  [38] 

E 

If  we  assume  values  likely  in  a  detection  scenario,  E  =  100  V/pm  (1  MV/cm),  d  =  15  pm  and  r  > 
10  ms,  the  minimum  mobility  required  is  therefore  5.0  x  10-7  cnr/V-s. 
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Figure  44.  Charge  collection  efficiency  and  detection  efficiency  as  a  function  of  detector 
thickness  for  ju  =  10-6  cnr/V-s,  E  =  1  MV/cm,  and  r  =  1  s. 


In  Figure  44,  charge  collection  efficiency  and  detection  efficiency  are  plotted  for  a  planar  thin- 
film  detector  with  /j,  =  10-6  cm/Vs,  x  =  1  s,  and  E  =  1  MV/cm.  The  charge  collection  efficiency 
increases  with  thickness  up  to  10  pm  (this  corresponds  to  the  distance  traveled  by  the  charge 
carriers  in  the  integration  time)  and  then  decreases.  The  initial  rise  is  due  to  an  increase  in  neutrons 
being  absorbed  near  the  incident  surface.  As  the  thickness  increases,  the  fixEld  ratio  decreases, 
which  leads  to  a  drop  in  charge  collection  efficiency.  The  detection  efficiency,  however,  continues 
to  increase  with  thickness,  proportionally  to  neutron  absorption  as  the  charge  collection  efficiency 
threshold  is  still  maintained  within  the  range  shown  in  the  graph. 
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Appendix  A 


Establishing  Electrical  Contacts  for  Amorphous  Hydrogenated  Boron  Carbide 

Michelle  M.  Paquette.1  Joseph  W.  Otto,1  Sudarshan  Karki,1  Bradley  J.  Nordell,1  M.  Sky  Driver,1  Thuong  D.  Nguyen,1  Christopher  L.  Keck,1  Brent  J.  Rogers,1 
Andrew  Kitahara,1  Sean  W.  King,2  William  A.  Lanford,3  Sudhaunshu  Purohit,4  Wenjing  Li,4  Chi  Zhang,4  Nathan  A.  Oyler,4  and  A.  N.  Caruso1 

Department  of  Physics  &  Astronomy,  University  of Missouri-Kansas  City,  Kansas  City,  MO,  United  States 
Logic  Technology  Development,  Intel  Corporation,  Hilsboro,  OR,  United  States 

Department  of  Physics.  University  at  Albany.  Albany.  NY.  United  States 
Department  of  Chemistry,  University  of  Missouri-Kansas  City,  Kansas  City,  MO,  United  States 


uivkc 

Department  ot  Physics 
and  Astronomy 


lu  limit  leakage 


(1)  To  maximize  device  efficiency  (e  g.,  ohnri 
or  recovery  of  charge,  or  blocking 

current/ minonry  carrier  injection) 

(2)  At  a  more  fundamental  level,  appropriate  electrical  contacts  are  critical 
to  studying  the  basic  electrical  ptopetties  of  a  thm  film  ot  device  (e.g.. 

Ohmic  contact  definition 

A  mctal-scmiconductof  device  with  ohmic  contacts  is  often  described  as 
having  a  linear  I  V  profile  fee.,  obeying  Ohm’s  law.  V  =  IR),  ot  as  having  low 
contact  resistance.  More  rigorously,  however,  we  can  say  that  “the  voltage 
drop  across  the  ohmic  onntacts  should  be  small  compared  to  the  voltage 
drop(s)  across  the  active  rcgion(s)  of  the  device,  or,  stared  another  way.  the 
ohmic  contact  should  supply  any  current  that  the  device  requires  in  its  normal 
mode  of  operation. “  (Shroder  et  at,  IEEE  Trans,  on  Electronic  Devices,  ED 
31,  1984,  637).  Importantly,  the  behavior  of  the  contact  cannot  he 
considered  in  isolation.  Whether  a  contact  may  be  defined  as  "ohmic"  or 
not  depends  on  how  it  behaves  relative  to  the  rest  of  the  device,  its  desired 

function,  and  the  operating  conditions. 

The  Schottky  model 

An  nveriy  simplistic,  yet  highly  useful  base  model.  Two  major  modes  of 

Transport  arc  thermionic  emission  (thermally  activated  carrier  transport  over 
the  barrier)  and  tunneling  through  the  barrier  in  the  case  ol  narrow  barriers. 


<o, 


bmuofimt  Lvrh.  Solitl  Surfaces,  Imfatcs  jM  Thm  Filets,  S'  El  Sprtuerf,  2C0t. 

Some  first-order  corrections  to  the  Schottky  model 

•  Application  of  bias  affects  bamcr  height. 

•  Carrier  concentration  (specifically  at  interface)  affects  barrier  width, 
interface, 

Some  corrections  to  the  Schottky  model  specific  to  amorphous  solids 

•  Interface  states  unlikely;  harrier  heighr  more  likely  to  correlate  with  metal 

•  Localized  states  within  the  gap  can  allow  hopping  transport  through  battier 

Considerations  specific  to  high- resistivity  solids 


limit  device  performance? 

•  Is  it  possible  to  form  perfect  blocking  CO 
localized  states  is  present  within  the  gap? 


taels  at  all  if  a  high  number  of 


Barrier  height  characterization 

•  Photocmission  spectroscopy 
(see  box  at  right). 

•  Can  use  in  conjunction  with 
electronic  structure  characterization 
methods  (e.g.,  UV/Vis  optical 

band  diagrams. 

•  I V  behavior  of  M-S  junction  (need 
Richardson's  constant) 

•  CV  profile  of  M  S  juncliun  (nol 
appropriate  for  high -resistivity 
materials) 

Electrical  behavior  of  contacts 
■  IV  ol'  MIM  structures 

(see  figure  at  far  right) 

•  Voltage  drop  across  contact 

(sec  figure  at  far  right) 

•  Contact  resistance 


Photoemission  Spectroscopy 


Inject 
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X-Ray  Photoemission  Spectroscopy 

Metal  overlayets  were  deposited  onto  a-B,C:H,  thin  films  to 
investigate  intctfactal  chemistry  and  electronic  structure  (see 
photocmission  ho*  above),  for  most  metals  snidics,  oxides  were 
observed  at  Ihe  interface,  even  though  the  metals  were  evaporated 

diffused  from  rhcbulk  a-B.C;H,  into  the  metals.  Example  results 
for  Pd  overlayers  ate  shown  below,  where  an  interfadal  layer  of 
rdO  is  fotmed  between  a-BrC:Il,  and  metallic  Pd.  Ko  band 
bending  was  apparent  in  this  case 
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Electrical  Characterization 

I  V  measuremenls  were  done  on  a  series  of  melal  insulator  metal  (MIM) 
hcteiosuuctures.  produced  under  both  in-situ  and  a-silu  conditions.  A  bottom 
contact  was  deposited  via  sputtering  onto  a  silicon  substrate,  followed  by  the 
deposition  of  an  a-B,C:H,  film  via  PECVD,  followed  by  deposition  of  a  top 
contact  via  sputtering.  For  ex-situ  conditions,  the  sample  was  exposed  to 
atmosphere  for  -30  min  between  depositing  one  or  both  of  the  layers. 


A  log-log  plot  of  one  branch  of 
the  above  data  displays  behavior 
typical  of  high-resistivity 
materials:  a  linear  region  at  low 
voltages  characteristic  of  ohmic 
behavior  and  a  linear  region  at 
higher  voltages  characteristic  of  a 
space-chargc-limired  current 


Samples  with  base  Cu  contacts  it 

handled  ex  slsu  display  more 
consistent  IV  behavior  across  ! 

voltage  ranges.  This  appears  to  be  £ 
a  reproducible  trend. 

Kv-situ  handling  docs  not  appear  o  ’ 


A  Cu/oxidc/a-B.CH/Pd  device 
displays  highet  current  and 
greater  linearity  than  the  all-Cu 


Hysteresis  curve  width  is  a 
function  of  step  size  (smaller  step 

size  =  narrower  curve). 


even  in  the  presence  of  an 
aluminum  oxide  layer  (very 
likely  Al,Ov  an  insulator) 
resulting  from  ex-situ  handling 
during  fabrication  Breakdown  in 
IV  behavior  occurs  at  >10  V, 
which  is  relatively  common 
among  the  devices  fabricated  to 


Summary  and  Conclusions 


■  Many  metals  are  expected  to  form  relatively  small  bartieis  at  the  a-B.C:H, 
interface;  initial  predictions  were  made  based  on  melal  work  functions. 

•  The  formation  of  imerfacial  oxides  due  ro  diffusion  of  oxygen  from  within 
the  bulk  a-B,C:H,  material  complicates  intcifacial  electronic  srruemre 
interpretation  and  may  create  additional  barriers  (or  decrease  barriers!) 

■  Cu  and  Td  appear  to  create  ohmic  contacts,  even  when  handled  under  a- 
situ  conditions,  on  the  basis  of  MIM  IV  curves.  However,  both  of  rhese 
metals  are  expected  to  form  accumulation-type  contacts  with  very  low 
barriers,  and  neither  oxide  is  expected  to  be  particularly  problematic. 

■  A1  also  appears  to  create  ohmic  contacts,  even  in  the  presence  of  an  oxide 
layer  (presumably  Al.Oj.  an  insulator).  This  suggests  that  the  higher  barrier 
and/or  insulating  oxide  aie  not  enough  to  create  significant  resistance  at 

the  contact  relative  to  the  rest  of  the  device,  and  therefore  the  contact  type 
may  not  he  particularly  important  for  fabricating  an  ohmic  device. 
However,  some  devices  definitely  exhibit  greater  integrity  at  >10  V,  which 
may  be  related  to  the  contacts. 

•  Although  the  analysis  of  MIM  IV  curves  may  be  useful  for  evaluating 
contact  resistance  relative  to  device  resistance,  it  may  be  too  simplistic  for 
evaluating  Ihe  quality  of  contacts  in  the  recovery  of  electron  hole  pairs  in 


Future  Work 


■  Continue  lo  explore  the  reproducibility  of  the  MIM  devices. 

•  Study  additional  metals  and  contact  materials  (Cr,  Ti.  Mi,  Pt,  ITO). 

•  Complete  more  rigorous  theoretical  analyses  on  the  MIM  I V  curves. 

•  Find  additional  methods  lot  investigating  contact  efficiency  in  the  context 
of  neutron  deteclois. 
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